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PREFACE 


The  collection  of  five  papers  included  in  this  report  are 
the  results  of  an  experimental  and  numerical  effort  undertaken  at 
Arisona  State  University  with  Office  of  Naval  Research  support 
under  contract  N00014-81-K-0428 .  The  paper  titles  and  where  they 
have  appeared  or  have  been  submitted  are  as  follows: 


1. 


U  Study  of  -Sudden  Expansion  -Pipe  -Plow  using  an 
Algebraic-Stress  Model  of  Turbulence*** by  B.K. 
Sultanian,  6.P.  Neitzel,  and  D.E.  Metzger,  AXAA  Paper 
No.  86-1062,  1986.  (Also  presented  at  the  1986 
AXAA/ASME  Fluid  Dynasties,  Plasma  Dynamics  and  Lasers 
Conference . ) 


2.  "Comment  on  'The  Flowfield  in  a  Suddenly  Enlarged 
Combustion  Chamber, A"  by  B.K.  Sultanian,  6.P.  Neitzel, 
and  D.E.  Metzger,  to  appear,  AXAA  Journal. 

3.  "Turbulent  Flow  Prediction  in  a  Sudden  Axisymmatric 
Expansion;"  by  B.K.  Sultanian,  6.P.  Neitzel,  and  D.E. 
Metzger,  in  Turbulence  Measurements  and  Flo*  Modmllng, 
Hemisphere  Publishing  Corporation,  New  York,  655-664, 
1987.  (Also  presented  at  the  1985  International 
Symposium  on  Refined  Flow  Modeling  and  Turbulence 
Measurements . ) 

4 .  "fieat  Transfer  to  Turbulent  Swirling  Flow  Through  a 
Sudden  Xxisymmetric  Expansion,,"  by  P.A.  Dellenback, 

D.E.  Metzger,  and  G.P.  Neitzel,  to  appear.  Journal  of 
Boat  Trans  far,  Trans »  ASMS.  (Also  to  be  presented  at 
the  1987  ASME/AXChE  National  Heat  Transfer  Conference.) 

i 

5.  -'"Measurements  in  Turbulent  Swirling  Flow  Through  an 
Abrupt  Axisymmetric  Expansion, "  by  P .A.  Dellenback, 

D.E.  Metzger,  and  G.P.  Neitzel,  submitted  to  AXAA 
Journal .  j 


The  first  three  papers  describe  numerical  predictions  of  the 
flowfields;/  paper- 4>  describes  experimental  measurements  of  local 
convective  haa^/t^rana f er  on  the  tube  walls  downstream  of 
expansion;  and/vpap4k:SX)  describes  detailed  mean  and  fluctuating 
flowfield  measurements  both  upstream  and  downstream  of  the 
expansion  including  processing  vortex  core  phenomena.  — 
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A  STUDY  OF  SUDDS*  EXPANSION  PIPS  PLOW  USING  AN  ALGEBRAIC 
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Indianapolis,  Indiana  46206 


6.  P.  Noitsol  and  D.  R.  Notz(ar 
Machanleal  and  Aorospaea  Enslnoarlns 
Arizona  Stats  Univarsity 
Tampa,  Arizona  85287 


Abstract 

The  pradictiva  ability  of  an  algebraic  stross 
nodal  (ASM)  and  tha  Boussinesq  viscosity  nodal 
(BVM),  or  tha  k-c  nodal,  is  tostad  on  s  suddan 
•spans ion  pips  flow  against  meant ly  published 
naasursnants  via  laser  Doppler  valociswtar  (LDV) . 
Calculations  ara  conparad  with  tha  naasursnants  on 
tha  naan  axial  velocity  and  axial  turbulence  inten¬ 
sity.  Whila  tha  standard  nodal  constants  are  used 
for  tha  BVM,  tha  additional  constants  that  nultiply 
the  "return- to- isotropy"  and  "rapid**  parts  of  tha 
pressure-strain  tern  in  tha  ASM  are  tuned  to  tha 
classical  pipe  flow  data  of  Laufer.  Both  In  tha 
recirculation  and  rodevelopnent  rations,  tha  ABM 
results  ara  in  excellent  at re an ant  with  tha  data, 
representing  an  inprovanent  over  tha  present  BVM 
sinulation  and  an  earlier  2/E/FIX  prediction  of 
this  flow.  Tha  ability  of  tha  ABM  to  successfully 
sinulate  tha  affects  of  streanllne  curvature  and 
anisotropy  in  tha  turbulence  field  appears  to  be 
tha  major  factor  contributing  to  this  success. 

Nonane  latum 


c  can tar line 

cl>  c2 
Ccl>  Cc2 

Cp  turbulence  nodal  constants 

H  step  height  •  R}  -  Rj 

in  inlet  _ 

k  turbulent  kinetic  energy  >0.5  ujuj 

out  outlet 

p  naan  pressure 

P  grid  point  (as  subscript) 

Q  volunetric  flow  rata 

Rj:  inside  pipe  radius  of  inlet  pipe 

*2  inside  pipe  radius  of  outlet  pipe 

S  source  tarn 

U£  naan- velocity  coaponent  in  x^  direction 

uj  fluctuating  velocity  coaponent  in  x^ 

dlmction 

u* ,v' ,w*  turbulence  intensity  coaponents 
x,y,r  coordinate  directions 

x^  coordinates  in  tensor  notation 

xr  reattaehaent  length 


•The  msearch  work  reported  hem  was  conducted  at 
Arizona  State  University. 
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diffusion  coefficient 

Krone cker  delta 

dissipation  rate  of  k 

ratio  of  production  to  dissipation  of  k 

molecular  viscosity 

eddy  (or  turbulent)  viscosity  -  cvpk‘/t 
klnanatic  eddy  viscosity  -  vt/p 
fluid  density 

turbulent  Prandtl  number  for  k 
turbulent  Prandtl  nuaber  for  c 
tenoral  dependent  variable 
belonging  to  variable  a  (as  subscript 
or  superscript) 


L _ mtrrtwUwi 

Engineering  consideration  of  a  sudden -expansion 
flow  geometry  dates  beck  to  1766  end  lorti'i* 
original  analysis.  However,  only  in  the  latter 
half  of  this  century  ties  the  problsa  received  a 
considerable  amount  of  interest.  The  aain  features 
of  a  pipe  flow  with  an  abrupt  expansion  am  de¬ 
picted  in  Figure  1,  which  shows  the  characteristics 
of  a  free-sheer  layer  near  the  expansion  and  those 
of  a  wall  boundary  layer  farther  downstream  of  the 
msttachmant  Near  the  expansion  face,  a  coherent, 
large-scale  turbulent  structure  plsys  an  iaportant 
role  in  the  subsequent  flow  development .  Although 
no  central  circulation  zone  exists  in  the  absence 
of  swirl,  the  existence  of  a  large  wall- bounded 
recirculation  region  makes  the  flow  predominantly 


Figure  1.  A  sudden-expansion  pipe  flow  structure. 
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elliptic  in  character.  Based  on  recent  LDV  met 
surements,  Stevenson2  confirms  the  presence  of  a 
counter- recirculating  eddy  within  a  half  step  H 
(difference  Rj-Rj  between  the  two  pipe  radii) 
of  the  expansion.  Note  that  at  the  beginning  of 
the  redevelopment  region  the  flow  has  many  features 
in  common  with  those  in  the  far  field  of  a  free- 
shear  layer,  quite  unlike  conditions  usually  found 
in  the  entrance  region  of  a  developing  pipe  flow. 
Lipstein3  reports  that,  in  the  range  of  B  •  0.4 
to  0.75  used  in  many  engineering  applications,  the 
reattachment  length  scales  properly  with  the  step 
height.  Perhaps  the  first  detailed  set  of  data  on 
both  the  mean  flow  and  the  turbulence  field  is  due 
to  Chaturvedi.4  However,  because  of  the  use  of 
an  intrusive  probe  (a  hot  wire)  for  measurement, 
the  data  are  not  completely  reliable  in  the  recir¬ 
culation  cone.  Using  these  and  other  experimental 
results,  Johnston3  presents  a  useful  discussion 
of  this  flow. 

A  flow- visualisation  study  by  Back  and 
Roschke,4  using  dye  injection,  demonstrates  the 
relative  Insensitivity  of  the  reattachment  length 
to  variations  in  Reynolds  number  in  the  fully  tur¬ 
bulent  flow  regime.  Drewry2  reports  an  Improved 
flow-visualisation  study  of  the  recirculation  re¬ 
gion  using  a  surface  oil-film  technique.  His  re¬ 
sults  show  that  the  reattachment  length  ranges  be¬ 
tween  7.9  and  9.2  step  heights  for  Reynolds  numbers 
(based  on  the  inlet  diameter)  in  the  range  of  1.3 
x  10*  to  2.2  x  10*.  The  studies  of  Hoon  and 
Rudinger8  and  Lu8  employ  LDV.  However,  their 
data  are  limited  to  the  mean  axial-velocity  distri¬ 
bution. 

As  to  the  past  computational  modeling  of  this 
elliptic  flow,  Cosman  at  al10  claim  to  accurately 
calculate  the  reattachment  length  observed  in  the 
experiments  of  Back  and  Roschke4  at  large  Rey¬ 
nolds  number.  Both  the  recirculation  region  and 
the  redevelopment  to  normal  pipe  flow  appear  to  be 
well  predicted  by  Ha  Hinh  and  Chassaing.33  Sayed 
and  Sturgess32  report  that  the  reattachment  length 
in  the  experiments  of  Chaturvedi4  is  underpradicted 
in  their  calculation  by  10-1SX.  Hoon  and  Rudinger8 
compare  numerical  results  with  their  LDV  data  for 
the  mean  flow;  good  agreement  between  the  two  is 
reported.  In  addition,  they  mention  that  changing 
one  of  the  model  constants,  i.e.,  Ctj  *  1.70  in- 
ataad  of  the  "standard"  value  of  1.92,  results  in 
better  agreement.  Most  of  the  results  just  men¬ 
tioned  have  been  obtained  using  the  TEACH  code  with 
the  standard  BVM.  The  flow  field  in  these  cases 
la  charactariced ,  at  the  inlet,  by  low  turbulence 
intensities  and  a  nearly  uniform  mean  axial-veloc- 
lty  profile,  with  the  exception  of  Moon  and 
Rudlnger8  where  a  fully  developed  inlet  velocity 
profile  is  used.  The  outflow  boundary  conditions 
in  all  these  cases  are  based  on  a  fully  developed 
(»/»x  ■  0)  flow  assumption. 

Host  of  the  gross  features  in  a  sudden  axisym- 
metrlc  expansion  are  also  shared  by  other  fully 
separated  internal  flows  such  as  those  over  single 
and  double  backward-facing  steps  and  in  confined 
jet  mixing  (ejectors).  An  overview  of  the  current 
capabilities  in  predicting  flows  in  this  class  are 
available  in  the  proceedings  of  the  1980-81  AFOSF- 
H1  In- Stanford  Conference.33  Sindir34  has  used 
an  ASM  as  one  of  the  four  models  for  predicting 
flows  in  backward- facing  step  geometries.  His 
study  shows  that  the  relative  performance  of  the 
models  is  region  dependent.  For  example,  the  "mod¬ 


ified"  ASH  has  been  found  to  produce  the-  best  pre 
dictions  in  the  reverse  flow  region  while  exhibit 
ing  too  slow  a  recovery  rate  beyond  the  reattach 
menl  tone. 

Recently,  Yang  and  Yu33  reported  a  complete 
set  of  LDV  data  for  an  isothermal  airflow  in  a 
dump- type  combustor  geometry.  The  flow  field  is 
somewhat  atypical  because  of  the  presence  of  a  high 
freestream  turbulence  level  at  the  inlet  and  a  con 
traction  at  the  outlet.  Computations  of  this  case 
using  BVM  found  that  the  mean  axial- velocity  data 
seriously  suffer  from  a  lack  of  mass  conservation 
downstream  of  the  expansion.  Professor  Yang33  sus¬ 
pects  that  the  error  could  be  due  to  a  shift  in 
the  blower  output.  A  strong  possibility  exists, 
however,  that  a  part  of  the  error  is  also  due  to 
velocity  bias,  which  is  known  to  worsen  in  regions 
of  high  turbulence  intensities.  The  importance  of 
such  a  bias  error  in  laser  velocimetry  is  clearly 
demonstrated  in  the  recent  investigation  of 
Stevenson  at  al.33 

The  measurements  of  mean  axial  velocity  and 
axial  turbulence  intensity  distribution  in  a  sud¬ 
den-expansion  pipe  flow  have  recently  been  reported 
by  Stevenson  et  al.38  These  are  very  accurate 
LDV  measurements  with  negligible  velocity  bias. 

The  present  calculations  of  both  the  mean  flow  and 
turbulence  field  based  on  the  turbulence  models 
BVM  and  ASM  are  compared  with  these  data  and  their 
aarlier  prediction,38  for  the  mean  axial  velocity 
distribution  only,  using  the  2/E/FIX  code  of  Pun 
and  Spalding.20 

L. _ Mathematical  Formulation 

LJ _ governing  Equations  with  ASM  and  BVM 

For  an  incompressible  turbulent  flow,  which  is 
statistically  stationary  in  the  siean,  the  conserva¬ 
tion  equations  for  mass  (continuity  equation)  and 
momentum  (Reynolds  equations)  may  be  written  as 
follows : 

8U. 


5^  ^jV  *  '  4  5^  (*  5^  '  'Vj)  (2> 


where  Uj  *  uj  is  the  local  instantaneous  velocity 
vector  and  -  p  ujuj  are  the  Reynolds  stresses. 
While  for  the  BVM,  the  Boussinesq  eddy  viscosity 
hypothesis , 


is  invoked,  for  the  ASM  (using  Rodi’s  hypothesis21) 
the  Reynolds  stresses  are  expressed  by  the  follow 
ing  algebraic  equations: 
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where  X  •  P^/c ,  the  ratio  of  local  production  of 
turbulent  kinetic  energy  to  its  diasipation.  The 
transport  aquations  for  k  and  c  are  given  by  the 
following: 


where 


#t  ■  vt^P  *  cv  k*/c ,  (7) 

and  ok  and  o(  are  the  equivalent  turbulent  Prandtl 
numbers  for  k  and  c,  respectively. 

Plow  predictions  are  greatly  influenced  by  the 
choice  of  constants  that  appear  under  any  turbu¬ 
lence  model  assumption,  lather  than  allowing  than 
to  be  used  arbitrarily  to  fit  data,  they  are  se¬ 
lected  with  hopes  of  having  universality.  The 
model  constants  used  in  the  BVM,  as  given  by 
Launder  and  Spalding,22  are  now  recognised  as 
standard  for  swat  flow  predictions.  Ro  such  gener¬ 
al  consensus  on  the  additional  modal  conatants  Cj 
and  C2  that  multiply  the  "return- to- isotropy"  term 
and  the  "rapid"  part,  respectively,  of  the  pres¬ 
sure-strain  tern  in  the  modeled  Reynolds  stress 
transport  equation23  has  yet  been  achieved,  al¬ 
though  Launder24  has  proposed  the  values  Cj  •  2.2 
and  C*  •  O.SS.  With  these  values,  the  ASM  re¬ 
sults*3  show  relatively  poor  comparison  with  the 
classical  pipe  flow  data  of  Laufer.2*  Since  pipe 
flow  is  the  asymptotic  flow  in  the  class  of  appli¬ 
cations  being  considered,  these  nodel  constants 
have  been  recalibrated23  to  yield  a  satisfactory 
pipe  flow  prediction.  Based  on  this  calibration 
study,  the  values  of  C}  •  2.2  and  Cy  -  0.70  have 
been  uaed.  Figure  2  depicts  the  comparison  of  mod¬ 
el  predictions  with  the  measurements  of  fully  de¬ 
veloped  velocity  profile  in  the  pipe. 

2.2  Boundary  Conditions 


Figure  2.  Comparison  of  calculations  with 
meaauremants  of  Laufer26  for  a  fully 
developed  pipe  flow. 


_ Solution  Procedure 


For  a  common  solution  procedure,  the  governing 
transport  equations  are  cast  into  the  following 
common  form: 


J  [{;  <•"■»>  •  b  >»«■»>]  • 

JfcKX)  *bKb)]*% 


(in 


For  a  general  elliptic  flow  in  a  pipe,  the  coatpu- 
tational  domain  has  three  types  of  boundaries: 
Inflow,  outflow,  and  no  flow  (solid  pipe  wall). 

The  flow  conaidered  here  permits  a  fully  developed 
(B/Bx  -  0)  boundary  condition  at  the  outlet.  Since 
the  BVM  and  ASM,  in  the  form  preaented  here,  are 
not  valid  in  the  region  of  low  turbulent  leynolds 
number  near  the  wall,  the  wall-function  approach 
of  Launder  and  Spalding22  la  uaed.  With  a  fully 
developed  flow  at  the  inlet,  while  the  one-seventh 
power-law  profile  ia  assumed  for  the  mean  axial 
velocity,  the  following  profilea  for  k*„  and  cln 
are  uaed: 

kln  •  0.0013  uj  (1  ♦  Afr/tj)2,5)  (8) 


where  q  represents  a  general  dependent  variable,  x 
and  r  art  the  axial  and  radial  coordinates  respec¬ 
tively  with  corresponding  velocity  components  U  and 
V,  r4  is  a  diffusion  coefficient,  and  is  the 
so-called  source  term.  The  latter  also  contains 
the  part  of  the  diffusion  term  that  cannot  be  ex¬ 
pressed  in  the  form  assumed  in  Equation  (11).  Fi¬ 
nite  difference  equations  for  each  q  are  obtained 
by  Integrating  Bquation  (11)  over  an  appropriate 
control  volume  constructed  sround  e  grid  point  rep¬ 
resenting  the  loestion  for  q.  The  details  on  de¬ 
riving  the  discretized  equations  in  this  way  are 
given  in  lef  27.  The  resulting  numerical  scheme  is 
based  on  the  "hybrid"  approximation  for  the  convec¬ 
tion  terms.  The  SIMPLE  procedure27  is  used  to  nu¬ 
merically  solve  these  equations. 


_3/4  .3/2„ 
*,  •  C  h .  /I 

in  v  in  1 


(9) 


where  the  mixing  length  1,  distribution  is  com¬ 
puted  from  the  equation 

Ig,  ■  1}  lO.lA-O.OKr/lj)2  -  0.0*(r/l1)4|  (10) 


The  convergence  problem  encountered  in  the  use 
of  ASM  has  been  handled  with  a  dual-loop  iteration 
scheme.23  According  to  this  scheam,  the  mean  flow 
and  turbulence  nodel  variables  (k  end  c)  are  solved 
in  an  outer  loop  using  Gauss-Seidel  line-by-line 
iteration  with  under-relaxation.  For  each  outer- 
loop  iteration,  the  leynolds  stresses  are  computed 
from  Equation  (4)  in  an  inner  loop  using  Gauss- 
Seidel  point  Iteration  with  under-relaxation. 


For  the  measurements  of  Stevenson.  Thompson, 
end  Could18  the  sudden  expansion  ipe  flow  geometry 
consists  of  pipes  of  diameters  76.2  mm  at  the  inlet 
and  152.4  mm  for  the  expansion  section  with  an  in¬ 
let  Reynolds  number  of  11  x  10*.  The  flow  is  fully 
developed  both  in  the  inlet  pipe  before  expansion 
and  at  the  exit  of  the  downstream  pipe  correspond¬ 
ing  to  x/H  «  40.  Although  both  the  biased  and  un¬ 
biased  LOV  measurement*  for  the  distribution  of 
swan  axial  velocity  and  turbulence  intensity  are 
reported  for  x/H  values  ranging  from  0.33  to  IS, 
only  the  unbiased  data  are  used  here  for  comparison 
with  the  computation.  In  the  experiment,  the  re- 
attaetment  is  found  to  occur  at  8.6  step  heights 
downstream  of  the  sudden  expansion. 

In  an  earlier  calculation19  of  this  flow  using 
the  2/E/FIX  computer  code  of  Pun  and  Spalding,20 
the  degree  of  agreement  with  data  worsens  as  one 
proceeds  downstream  in  the  redevelopment  region. 
Assuming  that  the  measured  velocity  profiles  are 
mass  conserving,  the  computed  ones  are  found  to  be 
deficient  in  this  respect.  In  this  calculation 
inlet  conditions  were  specified  at  x/H  «  0.33  (ig¬ 
noring  the  pipe  expansion  geometry  altogether)  and 
uniform  grid  spacing  was  used  in  the  axial  direc¬ 
tion,  which  might  not  properly  resolve  the  shear- 
layer  growth  in  the  recirculation  region.  With 
these  concerns  in  mind,  a  BVH  (with  its  standard 
model  constants)  prediction  using  a  SO  x  30  grid 
was  made.  The  grid,  shown  in  Figure  3,  is  nonuni- 
form  in  the  axial  direction  for  a  portion  of  the 
redevelopment  region  beyond  which  it  is  uniform  up 
to  the  outflow  boundary  at  x/H  •  40.  The  results 
for  the  mean  axial-velocity  distribution  are  com¬ 
pared  with  the  data  and  the  earlier  2/E/FIX  pre¬ 
diction  in  Figure  4.  Although  both  the  BVH  and 
2/E/FIX  predictions  themselves  are  in  good  agree¬ 
ment  in  the  recirculation  region,  excellent  agree¬ 
ment  between  the  data  and  BVH  prediction  is  seen 
at  x/H  •  15,  indicating  that  the  present  coaiputa- 
tion  is  mats  conserving  as  are  the  data. 


Figure  3.  Grid  system  50  x  30  used  for 
computing  the  flow  of  Stevenson  et  al.18 
(Grid  is  magnified  In  the  radial  direction.') 


good  agreement  with  the  data  in  the  recirculation 
region,  the  ASH  prediction  is  slightly  better  in 
the  redevelopment  region.  The  computational  re¬ 
sults  on  axial  turbulence  intensity  distribution 
for  the  two  models  are  summarised  in  Figure  6.  The 
profiles  have  their  maxima  along  the  extended  sur¬ 
face  of  the  inlet  pipe.  These  maxima  are  rather 
sharp  near  the  expansion  and  become  progressively 
flatter  downstream. 

From  the  results  for  both  the  mean  axial  veloc¬ 
ity  and  the  axial  turbulence  intensity  distribu¬ 
tions  it  is  seen  that  the  start  of  the  redevelop- 
ment  region  is  more  characteristic  of  the  far  field 
of  an  axisymmetric  jet  than  the  entrance  region  of 
a  conventional  pipe  flow.  This  is  evident  from  the 
fact  that  the  mean  axial  velocity  at  the  center 
line  decreases  from  x/H  -  9  to  x/H  -  15  and  the 
corresponding  turbulence  intensity  increases  over 
this  region.  In  a  conventional  pipe-flow  develop¬ 
ment,  turbulence  energy  is  produced  near  the  wall 
and  diffuses  toward  the  center;  whereas  in  a  sudden 
expansion  pipe  flow  just  the  opposite  occurs  in  the 
early  stages  of  redevelopment.  These  observations 
are  also  supported  by  the  experimental  investiga¬ 
tion  of  Chaturvedi*  and  its  discussion  by 
Johnston.8 


The  flow  calculation  with  the  ASH  alao  used  the 
grid  shown  in  Figure  3  with  boundary  conditions 
identical  to  the  ones  used  in  the  BVH  case.  The 
results,  alsr  shown  in  Figure  4,  are  in  excellent 
agreement  with  the  data  for  the  mean  axial  veloc¬ 
ity,  both  in  the  recirculation  and  redevelopment 
regions.  A  reattachment  length  of  8.5  atep  heights 
is  predicted  compared  with  the  experimentally  ob¬ 
served  value  of  8.6.  The  grid  has  an  inherent  er¬ 
ror*  of  2  x  10~3  and  about  500  iterations  are  re¬ 
quired  for  solution  convergence  according  to  the 
convergence  criterion28  used. 

In  the  present  calculations,  the  turbulence  in¬ 
tensities  are  equated  to  the  square  root  of  the 
corresponding  normal  Reynolds  stresses  divided  by 
the  maximum  mean  axial  velocity  at  the  inlet.  In 
the  case  of  ASH,  these  stresses  form  part  of  the 
solution  via  Equation  (4).  For  the  BVH,  however. 
Equation  (3)  is  used  to  compute  these  stresses. 
Results  for  the  axial  turbulence  Intensity  from 
both  the  BVH  and  ASH  predictions  are  compared  with 
the  deta  in  Figure  5.  While  both  results  are  in 

•Inherent  error  ■  I  (Qin  -  QoutW^in*  ''♦'•re  is 
the  known  flow  rate  at  the  inlet,  and  Qout  is  the 
computed  flow  rate  at  the  outlet  based  on  a  nomi¬ 
nal  fully  developed  velocity  profile  with  values 
at  the  grid  points. 


The  most  important  feature  of  an  internal  flow 
with  separation  is  the  recirculation  region  bounded 
by  the  dividing  streamline  and  the  confining  wall, 
as  shown  in  Figure  1.  At  the  recent  Stanford  con¬ 
ference,12  the  performances  of  several  turbulence 
models  in  computing  the  flow  in  a  backward- facing 
step  were  compared.  Host  of  the  BVH  versions  were 
found  to  underpredict  the  reattachment  length  by 
15-20%  when  compared  with  the  experiment  of  Kim, 
Kline,  and  Johnston.29  The  ASH,  however,  pre¬ 
dicts  reattachment  length  within  2%  of  the  experi¬ 
mental  value.  In  addition,  the  detailed  mean-flow 
results  in  this  case  are  in  excellent  agreement 
with  the  data,  especially  in  the  recirculation 
region. 

A  very  simple  model  is  suggested  in  Ref  30  for 
predicting  trends  in  the  reattaehment  length  data: 
"the  shear  layer  reattaches  when  it  has  entrained 
all  of  the  pressure  gradient  driven  backflow.” 

Since  the  pressure  gradient  is  primarily  dependent 
on  area  ratio,  accurate  prediction  of  reattachment 
length  depends  on  the  ability  of  the  turbulence 
model  to  simulate  the  entrainment  rate  (or  the 
growth  rate  of  the  shear  layer).  The  BVH  is  found 
deficient  in  this  respect,  predicting  a  higher 
growth  rate  for  the  shear  layer.  The  greatest 
weakness  of  all  second-order  turbulence  models 
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Figure  6.  Distribution  of  axial  turbulence 
intensity  computed  with  BVM  and  ASH  for  the 
flow  of  Stevenson  at  al.1® 


using  k-c  model  variables  is  the  transport  aqua 
tion  for  c  t  which  embraces  several  modeling  as¬ 
sumptions.  In  addition,  the  BVM  uses  a  simple 
eddy-viscosity  hypothesis.  To  partly  remedy  these 
drawbacks  and  to  capture  the  correct  growth  rate 
for  the  shear  layer,  the  sensitivity  of  reattach¬ 
ment  length  to  variations  in  the  model  constants 
Ct i  and  Ccj,  multiplying  the  production  and  dissi¬ 
pation  terms,  respectively  in  the  c-equation,  has 
been  studied.  Moon  and  Rudinger,®  for  example, 
recommend  changing  Ctj  from  its  standard  value  of 
1.92  to  1.70.  Stevenson  et  al,1®  based  on  their 
numerical  simulation  results,  report  a  linear  cor¬ 
relation  between  the  reattachment  length  and  C(2. 
Which  is  of  the  following  form: 

Cc2  «  0.059  <xr/H)  ♦  2.4305  (12) 


where  xr  is  the  reattschment  length.  Our  own  com¬ 
putation,  not  reported  here,  indicates  that  to  re¬ 
produce  the  experimentally  observed  reattachment 
length,  C(2  «  1.83  should  be  used  In  BVM.  Clearly, 
such  correlations  and  arbitrary  changes  in  model 
constants  have  little  general  validity  and  can  at 
most  be  justified  in  a  "postdictive"  computation. 

Cibson  et  al31  point  out  that  "the  turbulent 
shear  stress  and  intensity  are  reduced  by  stream¬ 
line  curvature  in  the  plane  of  the  mean  shear  when 
the  angular  momentum  of  the  flow  increases  in  the 
direction  of  the  radius  of  the  curvature  and  are 
increased  when  the  angular  momentum  decreases  with 
radius."  According  to  this  criterion  the  stream¬ 
lines  in  the  recirculation  region  of  a  pipe  expan¬ 


sion  flow  possess  a  stabilizing  curvature.  Even 
for  a  thin- shear  layer  prediction  the  utility  of  a 
linear  F- factor  correction  is  limited  to  flows  with 
mild  curvature.  So  and  Hellor3*-3*  >re  reasonably 
successful  in  simulating  the  effects  of  large  cur- 
vature  through  the  use  of  a  pressure- strain  term 
in  the  Reynolds  stress  equation.  Although  their 
method  is  open  to  criticism31  due  to  the  assumption 
of  local  equilibrium  and  the  omission  of  important 
mean-strain  effects  in  the  modeling  of  pressure 
strain  term,  the  approach  shares  the  general  con¬ 
clusion  drawn  by  Castro  and  Bradshaw3'  that  "meth 
ods  based  on  the  Reynolds- stress  transport  equa 
tions  will  be  needed  in  complex  flows."  The  con¬ 
traction  of  the  Reynolds  stress  transport  equations 
yields  the  equation  for  k;  the  pressure- strain  term 
becomes  identically  zero.  In  addition,  with  the 
eddy  -viscosity  assumption,  the  production  term  for 
k  involves  only  the  mean- strain  terms.  Thus,  the 
BVM  does  not  involve  any  term  representing  a  direct 
interaction  of  turbulence  with  the  mean  flow.  In 
this  respect  an  important  role  is  played  by  the 
"rapid"  part  in  the  ASM.  The  results  for  radial 
turbulence  intensity,  turbulent  kinetic  energy,  and 
turbulent  shear  stress,  shown  in  Figures  7  through 
9,  for  the  two  models  clearly  bring  out  the  superi¬ 
ority  of  ASM  over  BVM  in  successfully  predicting 
the  suppression  of  turbulence  by  stabilizing  curva 
ture  in  the  recirculation  region.  Further,  the 
strength  of  the  recirculation  eddy,  as  indicated 
by  the  maximum  value  for  the  stream  function  occur¬ 
ring  on  the  locus  of  flow  reversal,  predicted  by 
the  BVM  is  about  51  higher  than  that  predicted  by 
the  ASM. 
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Figure  7.  Distribution  of  radial  turbulence 
intensity  computed  with  BVM  and  ASM  for  the 
flow  of  Stevenson  et  al.1® 
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Figure  8.  Distribution  of  turbulent  kinetic 
energy  computed  with  BVM  and  ASH  for  the 
flow  of  Stevenson  at  al.18 


Figure  9.  Distribution  of  Reynolds  shear 
stress  computed  with  BVM  and  ASH  for  the 
flow  of  Stevenson  at  al.18 
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The  predictions  based  on  the  two  turbulence 
models  BVM  and  ASM  are  compared  with  a  recent  set 
of  LDV  measurements18  on  mean  axial  velocity  and 
axial  turbulence  intensity  in  a  sudden  expansion 
pipe  flow. 

The  ASH,  with  its  model  constants  tuned  to  a 
fully  developed  pipe  flow  data,  yielded  results  in 


excellent  agreement  with  the  measurements  in  both 
recirculation  and  redevelopment  regions.  The 
agreement  represents  an  improvement  over  the  BVM 
calculation  and  an  earlier  computation  of  this  flow 
using  the  2/E/FIX  code.  A  comparison  of  results 
for  the  turbulence  field  demonstrates  the  ability 
of  the  ASM  to  simulate  the  stabilizing  effects  of 
streamline  curvature  in  the  recirculation  region, 
which  results  in  relative  suppression  of  turbu¬ 
lence.  The  modeling  of  the  pressure-strain  redis¬ 
tribution  term  in  the  Reynolds- stress  transport 
equation  is  thought  to  be  accountable  for  this 
mechanism,  which  correctly  simulates  the  shear- 
layer  growth  and  hence  the  reattachment  length. 
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In  an  Interesting  paper,1  Professor  Yang  and  his  co-worker  Yu  report  a 
complete  set  of  laser  Doppler  veloclmeter  (LDV)  measurements  on  the  mean  veloc¬ 
ities  and  Reynolds  stresses  for  an  Isothermal  airflow  In  a  dump-type  combustor 
geometry.  A  maximum  error  of  6.5%  In  the  data.  In  terms  of  the  flow  rate  at 
each  measuring  station,  was  claimed  In  the  paper.  The  data  and  their  projected 

accuracy  motivated  us  to  use  them  to  validate  our  two-dimensional  elliptic  code 
2 

STEPUP  for  two  turbulence  models,  the  Bousslnesq  viscosity  model  (BVM),  or 
the  k-c  model,  and  the  algebraic  stress  transport  model.  Although  the  flow- 
field  Is  somewhat  atypical  because  of  the  presence  of  a  high  freestream  turbu¬ 
lence  level  at  the  Inlet  and  a  contraction  at  the  outlet,  computation  of  this 
flow  using  BVM  revealed  that  the  mean  axial-velocity  data  seriously  suffer  from 
a  lack  of  mass  conservation  downstream  of  the  expansion.  The  results  of  this 

♦This  Investigation  was  conducted  at  Arizona  State  University 
••Associate  Professor 
^Professor 


calculation  and  a  possible  explanation  for  the  error  In  the  data  are  briefly 
discussed  In  the  following  paragraphs. 

The  flow  geometry  In  this  case  corresponds  to  a  dump-type  combustor  where. 
In  view  of  a  short  chamber  length,  a  fully  developed  outflow  boundary  condition 
does  not  apply.  For  the  present  elliptic  computation,  specification  of  the 
measured  profiles  at  the  outlet  therefore  becomes  necessary.  Two  main  diffi¬ 
culties  are  encountered  In  trying  to  simulate  this  flow:  (1)  the  measured  mean 
axial-velocity  at  the  outflow  boundary  (x  «  40  cm)  yields  a  mass  flow  rate  that 
Is  about  65%  higher  than  that  at  the  Inlet,  and  (2)  the  turbulence  Intensity 

3 

at  the  Inlet  Is  unusually  high.  As  a  result,  the  commorriy  assumed  profiles 
at  the  Inflow  boundary  for  turbulent  kinetic  energy  k  and  Its  dissipation  rate 
c  are  found  to  be  Inadequate.  In  order  to  proceed  with  an  attempt  at  simula¬ 
tion,  the  mean  axial-velocity  profile  at  the  outlet  was  reduced  uniformly  to 
satisfy  continuity.  At  the  Inlet,  however,  along  with  a  nearly  uniform  mean 
axial  velocity,  the  following  profiles  for  k  (using  the  centerline  measurements 
of  turbulence  Intensities  reported  near  the  expansion)  are  used: 


kln  *  °*03  uln:  c1n 


except  at  the  near-wall  (Inlet  pipe)  grid  point  where 


(2) 


Is  specified.  The  form  In  Equation  (2)  Is  based  on  the  measurements  of 

4 

Holladay  In  developing  pipe  flow  that  show  that  the  maximum  value  of  k  near 
the  wall  Is  about  five  times  Its  value  In  the  free  stream  (core  flow).  The 


Inlet  condition  for  c. 


^In^P 


<k1n>P  '<0-41  V  • 


Is  based  on  the  existence  of  a  thin,  equilibrium  wall  boundary  layer. 


2 

Using  the  numerical  procedure  presented  elsewhere,  a  reattachment  length 
of  six  step  heights  (as  opposed  to  the  experimentally  observed  value  of  4.5) 
results  from  simulation  with  the  standard  BVM.  It  Is  seen  from  Figure  1  that 
the  computed  mean  axial -velocity  profiles  are  In  reasonable  agreement  with  the 
data  near  the  expansion  while  poor  quantitative  agreement  Is  obtained  further 
downstream.  Results  for  both  the  axial  and  radial  components  of  turbulence 
Intensity  along  the  chamber  axis  are  also  In  satisfactory  agreement  with  the 
measurements  shown  In  Figure  2.  Since  the  numerically  determined  mean  axial- 
velocity  profiles  are  mass  conserving,  the  disagreement  with  the  data  down¬ 
stream  of  the  expansion  Indicates  that  the  corresponding  experimental  profiles 
(In  addition  to  the  one  at  the  outflow  boundary)  are  also  Inconsistent  from  a 
continuity  consideration.  This  Is  further  verified  by  numerically  Integrating 
these  experimental  profiles  at  several  of  the  measurement  stations.  The  error 
at  each  station,  as  a  percentage  deviation  from  the  Inlet  flow  rate.  Is  also 
shown  In  Figure  2.  A  value  as  high  as  80%  for  this  error  can  be  seen  from 
Figure  2.  which  does  not  agree  with  the  reported1  value  of  only  6.5%. 

Professor  Yang5  suspects  that  the  error  could  be  due  to  a  shift  In  the 
blower  output.  A  strong  possibility  exists,  however,  that  a  part  of  the  error 
Is  also  due  to  velocity  bias  that  Is  known  to  worsen  In  regions  of  high  turbu¬ 
lence  Intensities.  The  Importance  of  such  bias  error  In  laser  veloclmetry  was 
clearly  demonstrated  In  the  recent  Investigation  of  Stevenson  et  al.5  In 
view  of  these  serious  discrepancies  In  the  data,  no  meaningful  validation  of  a 
numerical  procedure  nor  the  development  of  an  advanced  turbulence  model  are 
possible. 
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Figure  2.  Comparison  of  calculated  axial  and  radial  turbulence 
Intensities  with  measurements  at  the  pipe  centerline. 
Percentage  error  In  mass  flow  rate  computed  from  measured 
mean  axial-velocity  distribution. 
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SUMMARY 

Refined  modeling  of  sudden  expansion  pipe  flows  using  an  algebraic 
stress  transport  model  (ASM)  of  turbulence  is  presented.  The  model  has 
contributions  from  the  mean  flow  as  well  as  the  turbulence  field  In  Its 
pressure-strain  term,  and  the  corresponding  model  constants  are  tuned 
with  the  classical  fully-developed-plpe-flow  data  of  Laufer.  The 
results  of  computation  for  flows  with  Inlet  Reynolds  numbers  of  60,000 
and  110,000  are  compared  with  recent  L0V  measurements  and  also  with  the 
calculations  using  the  standard  Bousslnesq  viscosity  model  (BVM), 
alternately  known  as  the  k-c  model.  The  ASM,  with  its  ability  to 
naturally  simulate  the  effects  of  streamline  curvature  and  anisotropy 
In  the  turbulence  field.  Is  found  to  be  superior  to  the  BVM  In 
predicting  these  flows  In  both  recirculation  and  redevelopment  re¬ 
gions. 

1.  Introduction 

The  turbulent  flow  In  a  sudden  pipe  expansion  Is  an  Important  In¬ 
ternal  flow  with  separation  which  falls  In  the  general  class  of  complex 
shear  flows.  Such  a  flow  geometry  Is  of  common  occurrence  In  Industrial 
piping  systems  and  those  of  high-technology  aerospace  applications  with 
severe  volume  constraints.  A  dominant  feature  of  the  flow  Is  the  exis¬ 
tence  of  a  recirculation  zone  characterized  by  low  mean  velocities  but 
high  turbulence  Intensities.  A  typical  gas  turbine  combustor  flow 
field,  generally  three-dimensional  In  nature,  contains  one  or  more  such 
recirculation  regions  for  both  flame  stabilization  and  Increased  mixing 
of  fuel  and  air  In  Its  primary  zone. 

The  gross  features  found  In  a  sudden  axlsymmetrlc  expansion  are 
depicted  In  Fig.  1.  The  flow  structure  Is  rather  complex,  combining 
the  characteristics  of  a  free-shear  layer  near  the  expansion  and  those 
of  a  wall  boundary-layer  farther  downstream  of  the  reattachment.  The 
existence  of  a  large  wall-bounded  recirculation  region  makes  the  flow 
predominantly  elliptic  In  character.  The  early  part  of  the  redevelop¬ 
ment  region,  however.  Is  quite  different  from  the  usual  entrance  region 
In  a  developing  pipe  flow. 
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Fig.  1.  Flow  features  in  a  sudden  axisymmetric  expansion 


A  number  of  Investigations,  both  experimental  and  theoretical,  of 
a  sudden  expansion  pipe  flow  have  been  reviewed  In  Ref.  (14).  An  over¬ 
view  of  the  current  capabilities  In  predicting  flows  In  this  class  can 
be  had  from  the  proceedings  of  the  1980-81  AFOSR-HTTH-Stanford  Confer¬ 
ence  (3).  The  computational  work  reported  herein,  based  on  an  algebraic 
stress  model  (ASM)  and  the  standard  Bousslnesq  viscosity  model  (BVM) , 
uses  very  recent  LDV  data  In  this  flow  geometry  for  comparison. 
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For  an  Incompressible  turbulent  flow  which  Is  statistically  sta¬ 
tionary  In  the  mean,  the  conservation  equations  for  mass  (continuity 
equation)  and  momentum  (Reynolds  equations)  may  be  written  as 


where  Ui  +  ui  Is  the  local  Instantaneous  velocity  vector  and  -p  u-juj 
are  the  Reynolds  stresses.  The  ASM  used  here  Is  derived  from  the 
following  modeled  transport  equation  for  the  Reynolds  stresses 
presented  by  Launder  et  al.  (6): 
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The  pressure-strain  term  Is  split  Into  two  components  such  that 
(♦ij)-|t  which  represents  the  contribution  of  the  turbulence  field 
alone.  Is  modeled  using  the  "retum-to-lsotropy"  term  of  Rotta(ll)  and 
(+1j)2*  which  results  from  the  Interaction  of  the  mean  and  the 
fluctuating  field,  represents  the  "rapid"  part  whose  leading  term  Is 
retained  here. 

Invoking  Rodl's  hypothesis  (10)  results  In  the  following  algebraic 
equations  for  the  Reynolds  stresses: 


U1U1  -  I  k  *11  _ 


1  ",  C2  P1.1  "  3  Pk  *1j 


C-j  —  1  -*•  X 


(4) 


where  P|<  *  Pnn/2  an(1  x  *  P^/*,  the  ratio  of  local  production 
of  turbulent  kinetic  energy  to  Its  dissipation. 


For  BVM,  the  Bousslnesq  eddy  viscosity  hypothesis, 

_  aU1  •«!  2 

-  Vj  ■  ”t  +  57)  -  3  k 

Is  used.  The  transport  equations  for  k  and  c  are  given  by 
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and  tfk  and  ac  are  the  equivalent  turbulent  Prandtl  numbers  for  k  and  e, 
respectively. 


The  model  constants  used  In  the  BVM  are  those  given  by  Launder  and 
Spalding  (5)  which  are  now  recognized  as  standard  for  most  flow  pre- 
dlclons.  No  such  general  consensus  on  the  additional  model  constants 
Ci  and  C2  has  yet  been  achieved.  Presently,  the  values  of  -  2.2 
and  C?  *  0.70,  based  on  the  calibration  study  reported  In  (14)  against 
the  pipe  flow  data  of  Laufer  (4),  are  used. 
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For  a  general  elliptic  flow  In  a  pipe,  the  computational  domain  has 
three  types  of  boundaries:  Inflow,  outflow,  and  flow  (solid  pipe 
wall).  Both  the  flow  cases  considered  here  permit  a  fully-developed 
(3/ax  »  0)  boundary  condition  at  the  outlet.  Since  the  BVM  and  ASM,  In 
the  form  presented  here,  are  not  valid  In  the  region  of  low  turbulent 
Reynolds  number  near  the  wall,  the  wall-function  approach  of  Launder  and 
Spalding  (5)  Is  used  here.  With  a  fully-developed  flow  at  the  Inlet, 
while  the  one-seventh  power-law  profile  Is  assumed  for  the  mean  axial 
velocity,  the  following  profiles  for  k^n  and  c^n  are  used: 


kin  -  0.0013  U?  [1  ♦  4(r/R1)2-5j 


•In  -  Cl'* 


where  the  mixing  length  1„  distribution  Is  computed  from  the  equation 
(Schllchtlng,  12) 


lm  -  Rt  [0.14-0.08(r/R1)2  -  0.0B  (r/R^4) 


4.  Solution  Procedure 


For  a  common  solution  procedure,  the  governing  transport  equations 
are  cast  Into  a  common  form: 


r  tfe  <*"•♦>  *  Tr  '  r  &  <rr,  «>  *  f?  «rr,  £»  *  s,  <’*> 

»  A  »  a»n*r«1  rianandant  vaplahla  v  and  r  art  tha  avial 


where  +  represents  a  general  dependent  variable,  x  and  r  are  the  axial 
and  radial  coordinates,  respectively,  with  corresponding  velocity  compo¬ 
nents  U  and  V,  Is  a  diffusion  coefficient,  and  S^  Is  the 
so-called  source  term.  Finite  difference  equations  for  each  +  are 
obtained  by  Integrating  Eq.  (12)  over  an  appropriate  control  volume 
con-  structed  around  a  grid  point  representing  the  location  for  +. 
The  details  on  deriving  the  discretized  equations  In  this  way  are  given 
In  Ref.  (7).  The  SIMPLE  procedure  of  Patankar  and  Spalding  (8)  Is  used 
here  to  numerically  solve  these  equations. 


For  the  present  calculations  a  SO  x  30  (axial  x  radial)  non-uniform 
grid  Is  used.  The  convergence  problem  encountered  In  the  use  of  ASM  Is 
handled  with  a  dual-loop  Iteration  scheme  (14).  According  to  this 
scheme,  the  mean  flow  and  turbulence  model  variables  (k  and  c)  are 
solved  for  In  an  outer  loop  using  6auss-Se1de1  line-by-line  Iteration 
with  under-relaxation.  For  each  outer-loop  Iteration,  the  Reynolds 


4 


stresses  ere  computed  f row  Eq.  (4)  In  an  Inner  loop  using  Gauss-Sledel 
point  Iteration  with  under-relaxation.  The  outer-loop  Iterations  are 
terminated  according  to  the  convergence  criterion  based  on  the  'Inherent 
error*  of  the  computational  grid;  details  are  given  In  Ref.  (13). 

5.  Results  and  Discussion 

Parameter  values  for  the  two  sudden  expansion  pipe  flows  considered 
here  for  simulation  are  summed  up  In  Table  1.  For  Raj  «  60,000,  the 
calculations  of  the  mean  axial  velocity  distribution  with  the  two  turbu¬ 
lence  models  are  compared  with  the  L0V  data  In  Fig.  2  for  both  the  sepa¬ 
ration  (x/D]  »  1,  2,  3,  and  4)  and  redevelopment  (x/Oj  «  6,  8,  10, 
and  12)  regions  where  0^  Is  the  Inlet  pipe  diameter.  The  reattachment 
length  In  this  case  corresponds  to  x/D-|  -  4.5.  The  ASM-predlctlons 
are  seen  to  be  In  excellent  agreement  with  the  data  In  both  the  flow 
regions.  The  BMV-results,  on  the  other  hand,  are  significantly  differ¬ 
ent  from  the  experimental  data  over  the  latter  half  of  the  separation 
region  and  also  In  the  near-redevelopment  region. 


Table  1 .  Parameter  values  for  the  simulated  flows 


Flow 

R1 

*2 

Ul 

Re 

Data 

bTTB 

Fluid 

Hmrmi 

SHOalSH 

1 

source 

l 

Water 

■P9H 

49.3 

1.13 

60,000 

(1) 

2 

Air 

Bill 

76.2 

22.07 

110,000 

(2) 

For  the  second  case  with  Rei  -  110,000,  the  air  flow  Is  fully- 
developed  both  In  the  Inlet  pipe  before  expansion  and  at  the  exit  of 
the  downstream  pipe  corresponding  to  x/H  -  40  where  H  Is  the  step  height 
(R2  -  R^).  In  an  earlier  calculation  of  this  flow  by  Gould, 

Stevenson  and  Thompson  (2)  using  the  CHAMPION  2/E/FIX  computer  code  of 
Pun  and  Spalding  (9),  It  Is  found  that  the  degree  of  agreement  with  ex¬ 
perimental  data  worsens  as  one  proceeds  downstream  In  the  redevelopment 
region.  Since  the  measured  velocity  profiles  are  nearly  mass  conserv¬ 
ing,  the  computed  ones  are  found  to  be  deficient  In  this  respect  which 
has  been  attributed  to  a  possible  shortcoming  In  the  code  used. 

The  present  calculations  with  BVM  and  ASM  along  with  the  earlier 
2/E/FIX  prediction  for  the  mean  axial  velocity  are  compared  with  the 
data  In  Fig.  3.  It  Is  seen  that  although  both  the  BVM  and  2/E/FIX  pre¬ 
dictions  themselves  are  In  good  agreement  In  the  recirculation  region 
(x/H  ■  1,  3,  5,  and  7),  a  better  agreement  between  the  data  and  BVM 
prediction  Is  seen  In  the  far-redevelopment  region  (x/H  »  15,  Fig.  3). 
The  results  obtained  using  the  ASM,  however,  are  seen  to  be  In  excellent 
agreement  with  the  data  both  In  the  recirculation  and  redevelopment  re¬ 
gions.  With  this  model,  a  reattachment  length  of  8.5  step  heights  Is 
predicted  compared  to  the  experimentally  observed  value  of  8.6. 


i 


As  shown  In  Fig.  1,  the  most  Important  feature  of  the  flow  In  an 
sudden  axlsymmetrlc  expansion  Is  the  recirculation  region  bounded  by 
the  dividing  streamline  and  the  confining  wall.  An  accurate  prediction 
of  reattachment  length  depends,  among  other  things,  on  the  ability  of 
the  turbulence  model  to  simulate  the  entrainment  rate  (or  the  growth 
rate  of  the  shear  layer).  The  BVM,  with  Its  simple  eddy  viscosity  hypo¬ 
thesis,  Is  found  to  be  deficient  In  this  respect,  predicting  a  higher 
growth  rate  for  the  shear  layer. 

The  streamlines  In  the  recirculation  region  possess  a  stabilizing 
curvature  resulting  In  the  suppression  of  both  turbulent  shear  stress 
and  Intensity.  The  BVM,  using  the  model  variables  k  and  c  Is  unable 
to  accurately  simulate  this  vital  flow  physics.  This  Is  because,  the 
Reynolds  stress  transport  equations  upon  contraction  yield  the  equation 
for  k;  the  pressure-strain  term  becomes  Identically  zero.  In  addition, 
with  the  eddy-viscosity  assumption,  the  production  term  for  k  Involves 
only  the  mean-strain  terms.  Thus,  the  BVM  does  not  contain  any  term 
representing  a  direct  Interaction  of  turbulence  with  the  mean  flow.  In 
this  respect  an  Important  role  Is  played  by  the  "rapid"  part  (tij)? 

In  the  ASM.  The  results  for  turbulent  kinetic  energy  (Fig.  4)  and  tur¬ 
bulent  shear  stress  (Fig.  5)  for  the  two  models  clearly  bring  out  the 
superiority  of  ASM  over  BVM  In  successfully  predicting  the  suppression 
of  turbulence  by  the  stabilizing  curvature  In  recirculation  region. 

6.  Concluding  Remarks 

Two  sudden  expansion  pipe  flows  with  Inlet  Reynolds  numbers  of 
60,000  and  110,000  are  found  to  be  well  predlted  by  the  ASM  In  both 
their  recirculation  and  redevelopment  regions.  The  agreement  with  the 
data  for  these  flows  represents  a  significant  Improvement  over  the  pres¬ 
ent  BVM  calculations  and  an  earlier  computation  of  one  of  the  flows  us¬ 
ing  the  CHAMPION  2/E/FIX  code.  A  comparison  of  results  for  the  turbu¬ 
lence  field  using  both  the  BVM  and  ASM  demonstrates  the  ability  of  the 
latter  to  simulate  the  stabilizing  effects  of  streamline  curvature  In 
the  recirculation  region  which  results  In  relative  suppression  of  turbu¬ 
lence.  It  Is  believed  that  modeling  of  the  pressure-strain  redistribu¬ 
tion  term  In  the  Reynolds-stress  transport  equation  Is  accountable  for 
this  mechanism  which  correctly  simulates  the  shear-layer  growth  and 
hence  the  reattachment  length. 
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HEAT  TEARS FEB  TO  TUBBULENT  SWIBLIRG  FLOW 
THROUGH  A  SUDDEN  AXISTMNETBIC  EXPARSIOR 


P.A.  Dellenback,  D.E.  Metzger,  end  G.P.  Neitzel 
Arizona  State  University 


Experimental  data  are  presented  for  local  heat 
transfer  rates  in  the  tube  downstream  of  an  abrupt  2:1 
expansion .  Eater ,  with  a  nominal  inlet  Prandtl  number  of  6, 
was  used  as  the  working  fluid.  In  the  upstream  tube,  the 
Reynolds  number  was  varied  from  30,000  to  100,000  and  the 
swirl  number  was  varied  from  aero  to  1.2.  A  uniform  wall 
beat  flux  boundary  condition  was  employed,  which  resulted  in 
wall-to-bulk  fluid  temperatures  ranging  from  14*C  to  50*C. 
Plots  of  local  Russel t  numbers  show  a  sharply  peaked 
behavior  at  the  point  of  maximum  beat  transfer,  with 
increasing  swirl  greatly  exaggerating  the  peaking.  As  swirl 
increased  from  zero  to  its  maximum  value,  the  location  of 
peak  Russel t  numbers  was  observed  to  shift  from  8.0  to  1.5 
step  heights  downstream  of  the  expansion.  This  upstream 
movement  of  the  maximum  Russel t  number  was  accompanied  by  an 
increase  in  its  magnitude  from  3  to  9.5  times  larger  than 
fully  developed  tube  flow  values.  For  all  cases,  the 
location  of  maximum  heat  transfer  occurred  upstream  of  the 
flow  reattachment  point. 


HONE UCLA TORE 


Diameter  of  upstreaa  tube 
Diaaeter  of  downatreaa  tube 
Local  heat  transfer  coefficient 
Step  height;  (D2  -  D ) / 2 
Current  in  tube  wall 
Theraal  conductivity  of  water 
Mass  flow  rate 
Nuaselt  number 

Fully  developed  Nusselt  nuaber  for  turbulent 
flow  represented  by  Dittus-Boelter  or  Siede 
Tate  equations 

Maxiaua  or  peak  Nusselt  nuaber 
Local  heat  flux 
Total  heat  input 

Radial  position  relative  to  tube  centerline 
Resistance  of  tube  wall 

Reynolds  nuaber  in  upstreaa  tube;  Ud/v 
Upstreaa  tube  inside  radius 

Swirl  number  in  upstreaa  tube  as  defined  by 

2  U  V  dr 
U*  dr 

Bulk  fluid  temperature 
Outside  tube  wall  temperature 
Inside  tube  wall  temperature 
Local  mean  axial  velocity 

Axial  velocity  averaged  over  cross  section 
Maximum  axial  velocity  in  upstream  tube 
Local  mean  tangential  velocity 
Axial  distance  from  expansion  face 
Reattachment  length 

Axial  location  of  peak  Nusselt  no. 


turbulent  kinetic  energy  generated  there  dissipates 
relatively  slowly  maintaining  such  larger  levels  than  would 
be  found  in  ordinary  pipe  flow  where  no  such  internal  shear 
layer  exists.  With  high  levels  of  turbulence  kinetic 
energy,  diffusion  rates  are  elevated  and  the  thickness  of 
the  viscosity  dominated  sublayer  is  reduced,  resulting  in 
high  rates  of  heat  transfer  between  the  tube  wall  and  mean 
flow. 

Several  interesting  effects  appear  in  the  flowfield 
with  the  introduction  of  swirl  (Gupta  et  al.,  1984).  Among 
these  is  an  increase  in  growth  rate,  entrainment,  and  decay 
of  the  core  flow  just  downstream  of  the  expansion. 
Consequently,  the  flow  reattachment  zone  moves  upstream  as 
swirl  strength  is  increased.  Swirl  is  also  responsible  for 
increased  shear  rates,  greater  turbulence  production,  and 
longer  path  lengths  for  a  particular  fluid  particle  so  that 
the  effect  of  swirl,  like  the  effect  of  the  sudden 
expansion,  is  also  to  significantly  increase  heat  transfer 
rates  over  those  found  in  purely  axial  pipe  flow  (Hay  and 
West,  1975). 

A  further  complex  phenomenon  which  frequently  occurs 
in  swirling  flows  is  the  development  of  an  unsteady 
(although  usually  periodic)  asymmetry  in  the  flowfield. 
These  asymmetries  are  not  fully  understood,  but  it  is  known 
that  they  may  assume  many  different  forms  depending  on  flow 
geometry  and  swirl  strength  (Dellenback,  1986;  Leibovich, 
1984;  Hallett  and  Gunther,  1984).  At  low  swirl  levels  in 


the  present  study,  the  core  flow  departed  from  axial 
symmetry  in  the  neighborhood  of  the  expansion  and  then 
proceeded  to  precess  about  the  tube  centerline  (Dellenback) . 
This  flow  feature  is  referred  to  here  as  the  precessing 
vortex  core  (PVC)  after  Gupta  et  al.  (1984).  In  nost 
geometries  the  PVC  is  only  seen  for  swirl  strengths  large 
enough  to  produce  a  'bubble1  of  on-axis  recirculating  fluid, 
known  as  vortex  breakdown.  However,  in  the  present 
investigation  (and  in  the  work  of  Hallett  and  Gunther), 
oscillatory  flow  asymmetries  were  only  detected  at  swirl 
numbers  less  than  those  associated  with  vortex  breakdown. 
Specifically,  the  PVC  was  only  present  at  the  lowest  non¬ 
zero  swirl  number  for  each  of  the  three  Reynolds  numbers 
examined  in  this  study.  Flowfields  at  higher  swirl  numbers 
were  symmetric  with  no  apparent  PVC  (Dellenback). 

Present  computational  capabilities  are  only  able  to 
successfully  treat  some  of  the  simpler  limiting  cases  of  the 
present  problem  (Sultanian  et  al.,  1985).  Thus  the  purpose 
of  the  present  investigation  was  to  measure  heat  transfer 
characteristics  throughout  the  separation,  reattachment,  and 


redevelopment  regions  in  the  downstream  tube  for  a  variety 
of  flow  conditions. 

The  present  study  began  with  an  extensive 
characterization  of  the  turbulent  flowfield  throughout  the 
downstream  tube  and  at  two  locations  upstream  of  the 
expansion.  A  laser  Doppler  anemometer  (LDA)  was  used  to 
determine  axial  and  tangential  mean  velocities  in 


a)  Nu/Nufd  is  a  weak  function  of  Reynolds  number  for  air, 
hut  is  a  strong  function  of  Reynolds  nuaber  for  water. 

b)  Maxiaua  Nuaselt  nuabers  are  well  correlated  by 

Nua  =  C  Red2/3,  where  C  =  0.20  for  air  and  0.40  for 
water. 

c)  Locations  of  peak  Nuaselt  nuabers  aove  slightly  upstreaa 
with  increasing  expansion  ratio,  but  show  little 
dependence  on  Reynolds  nuaber. 

Recently,  Baughn  et  al.  (1984)  reported  an  extension 
to  Zeaanick  and  Dougall’s  (1970)  work.  In  previous 
experiaents,  the  region  just  downstreaa  of  the  expansion  has 
been  subject  to  relatively  high  rates  of  heat  conduction  in 
the  tube  wall  with  consequently  high  uncertainty  in  the  heat 
transfer  coefficients.  Baughn  et  al.  (1984)  devised  a  test 
section  specifically  to  ainiaize  the  axial  conduction 
prohlea,  thus  defeating  an  effect  which  they  suggest 
introduced  error  into  Zeaanick  and  Dougall’s  data.  For 
several  expansion  ratios,  a  ainiaua  in  the  Nusselt  nuaber 
was  observed  at  about  one  step  height  downstream  and  the 
authors  suggest  that  this  is  possible  evidence  for  a  very 
saall,  counter-rotating  corner-eddy. 

Habib  and  McEligot  (1982)  reported  an  ambitious 
calculation  of  the  flowfield  and  related  heat  transfer 
behavior  for  the  present  problem,  including  swirl.  Their 
published  results  appear  to  underpredict  the  peak  Nusselt 
nuaber  which  is  interpolated  from  the  data  of  Zemanick  and 
Dougall  (1970)  by  about  50*  for  an  unswirled  flow  with 


Re  =  50,000,  fi  =  0.5,  and  Pr  =  0.7.  However,  they  did  find 
the  reattachaent  zone  for  unswirled  flow  to  he  centered  at  a 
generally  accepted  value  of  8.2  atep  heights,  while  the 
corresponding  location  of  peak  Nusselt  number  occurred  at 
6.5  atep  heights. 


EXPERIMENTAL  APPARATUS 

A  stainless  steel  water  flow  loop  comprised  the  main 


element  of  the  test  facility. 


Swirl  was  generated  by 


tangential  slots  as  shown  in  Figure  1.  Inside  diameters  of 
the  axial  inlet  tube,  the  swirler  insert,  and  the  upstream 
test  section  were  5.08  cm.  The  sudden  expansion  was  12 
diameters  downstream  of  the  swirl  generator  and  the  axial 
inlet  tube  was  31  diameters  long  to  allow  axial  flow 
development.  Flowrates  to  the  slots  and  the  axial  inlet 
tube  could  be  controlled  independently,  thus  giving 
capability  for  a  continuously  variable  swirl  strength. 
Flowrates  of  the  tap  water  used  in  the  loop  were  measured 
with  turbine-type  flowmeters.  The  flow  loop  included  an  on¬ 
line  deaerator  for  removing  air  which  tended  to  leave 
solution  as  a  result  of  elevated  water  temperatures  near  the 
tube  wall. 

The  axial  and  tangential  components  of  mean  and  RMS 
velocities  were  measured,  one  at  a  time,  with  a  single- 
component  LDA.  The  profiles  were  taken  at  locations  both 
upstream  and  downstream  of  the  expansion  using  Plexiglas 


test  sections. 


Details  of  the  optical  system,  data 


analysis,  and  procedures  are  described  by  Dellenback  (1986) 


Following  the  flowfield  neasureaents,  heat  transfer 


tests  were  performed  in  a  horizontal  stainless  steel  tube  by 


passing  direct  current  in  the  tube  wall.  The  steel  heat 


transfer  test  section  was  of  the  same  dimensions  as  th< 


acrylic  test  section  used  for  the  velocity  measurements. 


Heating  began  uL,  ......  uui  include,  the  face  of  the 


sudden  expansion,  as  shown  in  Figure  2.  The  unheated 


upstream  tube  (inside  diameter  of  5.08  cm)  and  the  expansion 


face  were  the  same  Plexiglas  components  used  in  the 


flowfield  measurements. 


The  tube  downstream  of  the 


expansion  was  a  commercially  available  stainless  steel  tube 


with  an  inside  diameter  of  9.98  cm,  a  wall  thickness  of 


0.89  mm,  and  a  length  of  1.04  m.  Consequently,  the 


expansion  ratio  was  1.97:1.  Stainless  steel  flanges  were 


carefully  attached  to  each  end  of  the  tube  by  very  shallow 


welds  at  the  extreme  ends  of  the  tube.  Although  the  flange- 


to-tube  fit  was  snug,  it  was  assumed  that  all  of  the 


electrical  current  passed  through  the  weld  and  thus  power 


was  dissipated  over  the  full  1.04  m  length  of  the  tube. 


This  assumption  was  important  to  the  computation  of  heat 


flux  and  to  determining  the  location  of  the  point  of  maximum 


heat  transfer  since  the  flanges  were  relatively  thick  at 


1.27  cm  or  0.5  step  heights. 


Conduction  losses  from  the  heated  test  section  to  the 


upstream  tube  and  expansion  face  were  minimal  since  the 
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theraal  conductivity  of  Plexiglas  in  relatively  low.  Heat 
tranafer  to  the  nixing  plenua  iaaediately  downatreaa  of  the 
teat  aection  waa  niniaized  by  a  1.27  ca-thick  apacer  aade 
froa  a  aachinahle  dielectric  (Melaaine)  which  waa  placed 
between  the  teat  section's  flange  and  the  plenua.  Finally, 
nylon  bolta  were  uaed  to  aecure  the  teat  aection  to  the 
plenua  to  further  niniaize  conduction  heat  loaaea. 

Copper  electrical  huaaea  conaiating  of  an  inner  ring 
and  an  outer  ring  connected  by  six  spokes  were  aachined  froa 
a  single  plate  1.27  ca-thick.  The  busses  joined  current 
carrying  cables  to  the  stainless  steel  flanges  of  the  test 
section.  Experience  showed  that  a  liberally  applied  fila  of 
copper-based  anti-seize  coapound  between  the  bus  and  flange 
would  decrease  electrical  contact  resistance  and  consequent 
heat  generation.  A  dedicated  108  kW  direct  current  power 
supply  was  coupled  via  seven  2.5  ca-diaaeter  cables  to  the 
outer  rias  of  each  bus.  Two  water-carrying  copper  tubes 
were  soldered  around  the  periphery  of  the  busses’  outer  ria 
for  the  purpose  of  guard  heating  or  cooling.  Each  bus  had 
two  pairs  of  thermocouples  imbedded  in  such  a  way  that 
temperature  gradients  across  the  spokes  could  be  monitored. 
Due  to  the  large  power  levels  enployed  in  these  experiments 
(22  kW  at  8.5  V  and  2600  A),  ohmic  heating  in  the  cables 
proved  to  be  the  dominant  heat  source  in  the  power 
connection  scheme.  Hence,  guard  cooling  was  always 
performed  to  minimize  temperature  gradients  in  the  spokes  of 
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the  hue  and  heat  conduction  froa  the  cablea  to  the  bus,  end 
in  turn,  to  the  teat  aection. 

The  tube-wall  teiperature  diatributions  were  measured 
with  commercially  available  copper-constantan  thermocouples 
mounted  on  the  outside  of  the  tube.  The  Junction  of  each 
thermocouple  was  sandwiched  between  two  51  jta-thick  glass- 
reinforced  polymer-laminate  films  to  provide  electrical 


insulation  froa  the  test  section. 


A  high  thermal 


conductivity  paste  was  applied  between  the  thermocouples  and 
the  test  section.  Flat,  well  insulated  steel  bands 
encircling  the  test  section  held  the  thermocouples  tightly 
in  place.  Nineteen  thermocouples  were  employed,  spaced  at 
smaller  intervals  near  the  upstream  end  of  the  test  section 
to  provide  high  resolution  in  the  region  of  rapidly  changing 
heat  transfer  coefficients.  Four  inches  of  fiberglass 
insulation  surrounded  the  test  section  and  flanges,  and  all 
but  the  outer  rim  of  the  electrical  busses. 

Power  input  levels  to  the  working  fluid  were  deduced 
in  two  ways.  In  the  first,  power  was  taken  as  the  current- 
voltage  product  after  measuring  the  voltage  drop  across  both 
a  shunt  and  the  test  section.  A  thermocouple  was  affixed  to 
the  copper  shunt  to  correct  for  temperature  dependent 
resistivity.  Input  power  was  also  determined  from  the 
measured  fluid  enthalpy  rise  through  the  test  section.  To 
this  end,  bulk  fluid  temperatures  were  measured  at  the  inlet 
and  outlet  of  the  test  section  with  immersion  thermocouples. 


iV. 


Over  the  range  of  test  conditions,  wall-to-bulk  fluid 
temperatures  ranged  from  14* C  to  50* C. 

PROCEDURES  AND  DATA  REDUCTION 

The  22  kW  power  levels  produced  heat  fluxes  on  the 
order  of  6.7  W/cm2 .  After  applying  power  to  the  test 
section,  the  time  required  for  the  loop  to  reach  steady 
state  was  between  30  and  60  minutes. 

For  most  trials,  the  wal 1- temperature  measuring 
thermocouples  were  located  along  the  top  of  the  tube. 
However,  in  order  to  examine  the  influence  of  natural 
convection  on  the  present  problem,  for  at  least  one  of  the 
multiple  trials  at  each  flow  condition  the  entire  test 
section  (with  thermocouples  still  attached)  was  removed, 
rotated  180*  about  the  tube  axis,  and  reinstalled  so  that 
the  temperature  measurements  could  be  made  along  the  bottom 
of  the  tube. 

Local  heat  transfer  results  are  presented  in  terms  of 
Nusselt  numbers  normalized  with  those  for  fully  developed 
non-swirling  flow,  where 


Here,  D2  is  the  downstream  tube  diameter,  k  is  the  thermal 
conductivity  of  water  evaluated  at  the  local  film 
temperature,  and  h  is  the  local  heat  transfer  coefficient 


The  local  heat  flux  is  designated  by  q,  the  local  wall 
tesperature  by  TM ,  and  the  local  bulk  teaperature  by  Tb . 
Because  wall-to-bulk  fluid  teaperature  differences  were 
moderate-to-large,  and  hence  property  variations 
appreciable,  the  Sieder-Tate  correlation  (Kern,  1950)  was 
used  to  evaluate  reference  Nusselt  numbers  for  the 
corresponding  fully-developed  flow  without  swirl: 

Nufd  =  0.027  He  •  Pr*'2  (3) 
All  fluid  properties  used  in  eqn.  (3)  were  evaluated  at  the 
local  bulk  teaperature  except  (t„  ,  which  was  evaluated  at  the 
local  wall  teaperature. 

Although  the  present  experiaents  have  a  nominally 
unifora  heat  flux  boundary  condition,  the  local  heat  fluxes 
are  not  strictly  unifora  due  to  both  the  teaperature 
dependence  of  tube  aaterial  properties  and  axial  heat 
conduction  in  the  tube  wall.  The  procedure  used  for  finding 
the  local  heat  fluxes  began  by  dividing  the  wall  into  19 
control  volumes  corresponding  to  the  19  wal 1- temperature 
theraocouples .  The  magnitude  of  Joule  heating  in  each 
volume  was  determined  from  I2R,  where  a  linear  curve  fit  of 
the  teaperature  dependent  resistivity  of  the  stainless  steel 
(Touloukian,  1967)  test  section  was  used.  Heat  loss  from 
the  outside  of  the  tube  wall  through  the  fiberglass 
insulation  which  surrounded  the  test  section  was  assumed 
negligible.  This  loss  was  eventually  calculated  to  be  no 
aore  than  0.05%  of  the  total  power  input  in  a  worst  case 


analysis  of  the  data.  The  wall  teaperature  (Tw)  in  eqn.  (2) 
is  the  teaperature  of  the  inside  surface  of  the  teat  section 
wall.  Because  exterior  surface  teaperaturea  were  actually 
measured,  the  interior  surface  teaperaturea  were  deterained 
froa  the  solution  of  the  one-diaensional  heat  conduction 
equation  for  a  cylindrical  shell.  Thia  procedure  aasuaea 
that  only  radial  conduction  is  important,  an  assuaption  that 
auat  be  justified  since  natural  convection  nay  cause  slight 
asyaaetry  in  circuaferential  teaperature  profiles.  Also, 
the  axial  variation  in  heat  transfer  coefficients  gives  rise 
to  axial  teaperature  gradients  and  conduction  in  the  tuhe 
wall.  However,  both  circuaferential  (discussed  below)  and 
axial  teaperature  gradients  were  exaained  and  found  to  be 
saall  coapared  to  those  in  the  radial  direction.  A  worst 
case  coaputation  of  axial  heat  flux  for  the  present  results 
showed  that,  at  aoat,  the  net  heat  conducted  out  of  a  tube 
wall  elenent  to  an  adjacent  element  was  about  IX  of  the 
Joule  heat  generated  in  that  element,  even  near  the  test 
section  flanges. 

Determining  the  local  heat  flux  and  inside  wall 
temperature  is  inherently  an  iterative  process,  but  one 
which  converges  suitably  in  one  iteration  because  the 
teaperature  differences  across  the  tube  wall  are  small 
(typically  on  the  order  of  2*C).  The  maximum  variation  in 
local  heat  flux  between  any  two  locations  od  the  tube  was 
eventually  calculated  to  be  less  than  1.5X.  Finally,  to  he 
consistent  in  methodology,  local  water  enthalpies  and  bulk 


a? 


temperatures  were  computed  from  an  energy  balance  which 
considered  the  local,  slightly  non-uniform,  heat  input. 

The  methods  of  Kline  and  McClintock  (1953)  were 
employed  to  determine  that  the  largest  uncertainties  are 
about  2%  in  upstream  Reynolds  number  and  8X  in  swirl  number. 
The  highest  uncertainties  in  Nusselt  numbers  were  computed 
to  be  about  9%,  with  these  occurring  near  the  location  of 
peak  Nusselt  number  where  wall-to-bulk  temperature 
differences  were  smallest.  Details  of  the  analysis  and 
tabular  data  are  given  in  the  dissertation  by  Dellenback 
(1986) . 

RESULTS  AND  DISCUSSION 

In  the  figures  which  follow,  the  axial  coordinate  is 
non-dimensional ized  with  the  upstream  tube  diameter  (D). 
The  swirl  numbers  shown  on  the  plots  were  evaluated  from 
velocity  profiles  measured  two  diameters  upstream  of  the 
expansion.  Similarly,  Reynolds  numbers  are  based  on 
properties  and  mean  velocities  in  the  upstream  tube. 

Summary  of  the  Flowfield.  Figures  3  and  4  show 
typical  mean  velocity  profiles  for  various  swirl  levels  at 
Re  =  100,000.  All  mean  velocities  (and  the  turbulence 
intenaities  referred  to  below)  were  normalized  using  the 
maximum  axial  velocity  in  the  upstream  tube.  Data  for 
unswirled  flows  agreed  very  well  with  the  recent  LDA  data  of 
Stevenson  et  al.  (1983),  but  flow  predictions  of  Habib  and 


McEligot  (1982)  for  strongly  swirled  flows  are  in  only  fair 


agreement  as  to  locations  of  peak  velocities  near  the 


expansion  and  extents  of  recirculation  regions. 


Flow  conditions  upstreas  of  the  expansion  varied 


greatly  with  swirl  nuaber.  Upstream  velocity  profiles  for 


the  lowest  swirl  case  at  each  Reynolds  nuaber  closely 


approxiaated  solid  body  rotation  plug-flow.  Figure  3  shows 


that  with  increasing  swirl  nuaber  the  highest  axial 


velocities  aove  toward  the  tube  wall.  Near  the  tube 


centerline,  large  gradients  in  the  aean  tangential  velocity 


are  responsible  for  turbulence  intensities  on  the  order  of 


30*.  Turbulence  diffuses  away  froa  the  tube's  axis  with 


both  the  axial  and  tangential  turbulence  intensities 


decreasing  rapidly  to  about  12*  between  r/R  =0.2  and  the 


wall . 


Prandtl  numbers  of  the  fluid  entering  the  heated  test 


section  were  5.1  at  Re  =  100,000;  5.8  at  Re  =  60,000;  and 


from  5.9  to  6.7  for  Re  =  30,000.  This  variation  in  Prandtl 


number  is  probably  insignificant  when  examining  the  present 


results,  as  suggested  by  the  near  congruence  of  Krall  and 


Sparrow's  (1966)  data  for  Prandtl  numbers  of  3  and  6. 


Axial  Heat  Transfer  Variations.  Figures  5-7  show  the 


measured  axial  variations  in  normalized  Nusselt  number  as  a 


function  of  swirl  for  nominal  Reynolds  numbers  of  30,000, 


60,000,  and  100,000,  respectively.  The  peak  Nusselt  numbers 


increase  consistently  in  magnitude  and  move  upstream  with 


increasing  swirl  strength.  This  upstream  migration  of  Nua 


is  a  direct  result  of  the  shortening  reattachment  length 
discussed  previously.  The  shortening  of  the  recirculation 
region  causes  shear  rates  and  hence  production  of  turbulence 
kinetic  energy  to  increase  with  consequently  higher  heat 
transfer  rates.  This  enhanceaent  is  also  proaoted  by  higber 
local  aean  velocities  as  the  tangential  velocity  coaponent 
increases  at  noainany  uuunaul  values  of  aean  axial 
velocity. 

Coaparison  of  the  unswirled  flow  results  fron 
Figures  5-7  deaonstrate  that  larger  enhanceaents  in  heat 
transfer  rates  over  straight  pipe  flow  occur  at  lower 
Reynolds  nuabers.  This  feature  is  siaply  rationalized  by 
the  observation  that  convection  heat  transfer  behavior  in 
separated  flow  is  coaaonly  found  to  depend  on  Re2/3  (Krall 
and  Sparrow,  1966;  Zeaanick  and  Dougall,  1970),  so  that  when 
Nusselt  numbers  are  normalized  with  a  fully-developed  tube- 
flow  value  which  depends  on  Re0-8  ,  the  ratio  depends 
on  Re-013.  For  zero  swirl,  the  recovery  to  fully-developed 
flow  Nusselt  numbers  occurs  faster  with  increasing  Reynolds 
number.  Nusselt  numbers  are  strongly  influenced  by 
turbulence  intensity,  which  velocity  field  measurements  gave 
as  a*  BX  for  Re  =  30,000  and  as  3k  for  Re  =  100,000  at  the  end 
of  the  test  section.  The  turbulence  intensities  have 
relaxed  more  quickly  in  the  latter  case  due  to  higher  rates 
of  dissipation  at  the  larger  Reynolds  numbers  (Laufer, 
1954) . 


Table  1  gives  a  auiiary  of  peak  Nuaaelt  numbers  (Nu,), 
locations  of  NuB  (x/H„u  ),  and  reattachment  lengths  (xr/H) 
for  the  flows  examined  in  this  study.  Because  there  was 
an  inevitable  limit  to  the  spatial  resolution  of  the 
temperatures  which  could  be  attained,  magnitudes  and 
locations  of  the  maximum  Nusselt  numbers  which  are  tabulated 
in  Table  1  were  determined  by  interpolation  of  expanded 
scale  plots  of  local  Nusselt  numbers  in  this  region.  The 
tabular  results  show  that  the  locations  of  maximum  Nusselt 
number  are  largely  Reynolds  number  independent,  consistent 
with  the  unswirled  results  of  investigations  discussed 
previously. 

Maximum  Nusselt  numbers  are  shown  for  the  present  data 
and  five  other  investigations  in  Figure  8.  The  present 
unswirled  data  for  NuB  fit  comfortably  among  similar  water 
data  which  show  the  Re2/3  dependence  mentioned  above. 
However,  at  the  larger  swirl  numbers  there  is  a  slight,  but 
noticeable  decrease  in  the  slope  of  the  curves.  Although 
the  swirl  number  is  not  a  fixed  value  at  each  Reynolds 
number,  we  may  tentatively  infer  from  Figure  8  that  NuB 
depends  approximately  on  Re0-46  for  a  fixed  swirl  number  in 
the  range,  0.60  <  S  <  1.23.  Thus  it  appears  that  heat 
transfer  enhancement  due  to  swirl  alone  does  contain  some 
Reynolds  number  dependence. 

Figure  9  contrasts  the  locations  of  maximum  Nusselt 
number  with  the  flow  reattachment  points  determined  from  LDA 
measurements.  Although  it  is  often  assumed  that  the  peak 
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Nusselt  number  occurs  at  the  reattachment  point  (Krall  and 
Sparrow,  1966;  Zemanick  and  Dougall,  1970),  for  this  work  it 
fell  consistently  upstream  of  the  reattachment  point  for  all 
swirl  numbers.  The  axial  distance  is  non-dimensionalized 
with  the  step  height  (H)  of  the  expansion  in  Figure  9  since 
reattachment  lengths  have  been  frequently  shown  to  correlate 
well  with  this  length  scale. 

Maximum  Nusselt  numbers  display  a  consistent  behavior 
as  swirl  number  is  varied  at  constant  Reynolds  number,  as 
shown  in  Figure  10.  These  curves  lead  to  speculation  that 
swirl  influence  on  heat  transfer  may  be  expressible  as  a 
power  law,  as  for  example  in  the  work  of  Hay  and  West 
(1975),  where  Nu/Nufd  =  (S  +  l)1-75  for  free  swirling  flow 
in  a  constant  diameter  pipe.  However,  an  indirect  influence 
of  the  sudden  expansion  is  that  it  causes  the  swirl  number 
to  change  dramatically  Just  downstream  of  the  expansion 
(Gupta  et  al . ,  1984;  Dellenback,  1986),  and  this  feature, 
along  with  the  diff icult-to-specify  relationship  between  Nua 
and  Reynolds  number,  have  thus  far  defeated  efforts  to 
separate  quantitatively  the  effects  of  separation  and  swirl 
in  the  present  results. 

Computed  differences  between  the  two  methods  which  are 
available  for  determining  energy  balances  are  also  shown  in 
Table  1  (as  *X  difference  in  Q’).  Incomplete  mixing  in  the 
downstream  plenum  is  believed  to  be  largely  responsible  for 
the  discrepancies  by  giving  larger  than  actual  downstream 
bulk  temperatures.  This  hypothesis  is  consistent  with  the 


mi 


Jinx’, 


downstream  thermocouple's  location  which  would  be  within  the 


wall  boundary  layer  if  the  test  section  extended 


approximately  30  cm  into  the  plenum.  In  almost  every  case, 


Z(I2R)  was  smaller  than  mCpATb,  which  is  also  supportive  of 


the  present  hypothesis  but  contrary  to  normal  expectation, 


where  heat  losses  usually  cause  the  fluid  enthalpy  rise  to 


he  less  than  the  power  input. 


Rapidly  fluctuating  wall  temperatures  were  observed  at 


low  Reynolds  and  swirl  number  combinations,  indicating  a 


random  unsteadiness  at  these  flow  conditions.  For  example, 


at  Re  =  30,000  (for  both  S  =  0  and  0.14)  wall  temperatures 


were  repeatedly  observed  to  fluctuate  by  as  much  as  14*C 


over  a  short  time  interval  of  about  two  seconds,  but  these 


fluctuations  were  not  periodic.  Similarly,  for  Re  =  60,000 


(both  S  =  0  and  0.18),  fluctuations  were  on  the  order  of 


*2*C  for  the  same  time  frame. 


When  these  fluctuations 


occurred  at  the  lowest  swirl  strength,  there  was  no  apparent 


correlation  with  the  PVC  previously  mentioned.  Also, 


because  the  unsteadiness  was  apparent  for  the  unswirled  flow 


as  well  as  the  weakly  swirled  flow,  it  is  plausible  that  the 


fluctuations  are  a  consequence  of  slight  shifts  of  the 


reattachment  point  within  a  small,  but  finite,  zone  of 


reattachment.  At  higher  Reynolds  and  swirl  numbers  (for 


which  the  reattachment  zone  i9  narrower),  these  large 


fluctuations  disappeared. 


Impact  of  Free  Convection.  The  influence  of  free 


convection  on  heat  transfer  in  these  experiments  was 


contribution  by  free  convection  occurred  at  the  last  four 


downstream  measurement  stations  for  Re  =  30,000.  For 
consistency  in  comparing  low  swirl  and  low  Reynolds  number 
cases,  only  runs  where  all  temperature  measurements  were 
made  using  the  bottom-mounted  thermocouples  have  been  used 
in  this  presentation. 


Comparison  with  Previous  Investigations.  Comparison 
of  the  present  results  with  those  of  other  investigations  is 
limited  to  the  case  of  unswirled  flows  for  which  data  are 
readily  available.  Figure  11  shows  one  such  comparison  for 
86,400  <  Re  <  101,000.  The  comparison  is  based  on 
normalized  Nusselt  numbers,  hut  since  fully-developed 
Nusselt  numbers  used  for  normalization  were  obtained  in  a 
different  manner  for  each  study,  this  can  present  problems 
when  comparing  results.  For  example,  if  the  present  results 
had  been  normalized  using  the  0 i t tus-Boel ter  relation  rather 
than  the  Sieder-Tate  correlation,  then  peak  values  of 
Nu/Nufd  would  have  been  8k  larger  for  Re  =  100,000  (and 
about  Ilk  larger  for  Re  =  30,000)  than  the  values  shown  in 
Figures  5-7  and  11.  However,  a  study  by  Malina  and  Sparrow 
(1964)  concluded  that  the  Sieder-Tate  relation  over-predicts 
variable-property  enhancement  of  heat  transfer  by  about  4k 
while  the  Dittus-Boelter  relation  undei — predicts  fully- 
developed  tube  flow  Nusselt  numbers  by  6-llk  in  the  Reynolds 
and  Prandtl  number  ranges  of  the  present  data.  Thus  the 
results  of  Malina  and  Sparrow  suggest  that  the  apparent  10k 
difference  in  Nu/Nuf(J,  which  is  due  to  the  choice  of  a 
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correlation  for  normalization,  nay  in  fact  be  much  aaaller. 
For  the  work  of  Krall  and  Sparrow  (1966),  the  data  of  Malina 
and  Sparrow  were  used  to  obtain  fully-developed  Nuaselt 
numbers,  but  the  effect  of  eaploying  the  Dittus-Boelter 
relation  would  be  to  increase  the  peak  Nu/Nuf(J  ratio  by  5- 
10*  over  their  values  shown  in  Figure  11.  To  complicate  the 
comparison  further,  Zemanick  and  Dougall  (1970)  used  their 
own  downstream  data  for  normalization,  while  the  data  of 
Baughn  et  al .  (1984)  and  Ede  et  al.  (1956,  1962)  are 
normalized  with  the  Dittus-Boelter  relation.  In  any  event, 
the  Nusselt  numbers  of  Krall  and  Sparrow’s  appear  somewhat 
high  at  several  Reynolds  numbers.  It  seems  reasonable  to 
conclude  that  the  vena-contracta  produced  by  their  orifice 
makes  the  effective  expansion  ratio  larger  than  2:1,  and  it 
has  been  clearly  demonstrated  (Krall  and  Sparrow;  Zemanick 
and  Dougall;  Baughn  et  al.)  that  peak  Nusselt  numbers 
increase  with  increasing  expansion  ratio.  The  present 
results  exhibit  very  good  agreement  with  the  recent  results 
of  Baughn  et  al.,  particularly  at  small  x/D.  This  provides 
confidence  in  the  present  results  since  the  principle 
motivation  for  the  work  of  Baughn  et  al.  was  to  minimize 
axial  heat  conduction,  especially  in  the  near  expansion 
region  where  measurement  inaccuracies  are  usually  largest. 


CLOSDBE 


The  principal  contribution  of  the  present  results  is 
to  specify  the  quantitative  effect  of  swirl  on  the  heat 
transfer  enhancement  for  the  flow  downstreas  of  a  sudden 
axisymmetric  expansion.  Local  heat  transfer  rates  have  been 
shown  to  increase  dramatically  in  the  separated  flow  region 
and  to  peak  just  upstream  of  the  flow  reattachment  point. 
The  results  also  show  that  the  location  of  peak  Nusselt 
number  is  a  strong  function  of  swirl  number.  However,  the 
locations  of  maximum  heat  transfer  rate  are  found  to  be 
largely  Reynolds  number  independent,  particularly  at  higher 
swirl  numbers. 
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MEASUREMENTS  IN  TURBULENT  SWIRLING  PLOW 
THROUGH  AN  ABRUPT  AXISYMMBTRIC  EXPANSION 

P.A.  Dellenback,  O.E.  Metzger,  and  G.P.  Neitzel 
Arizona  State  Univeraity 


Experimental  data  are  presented  for  both  axial  and 
tangential  velocity  components  in  turbulent  swirling  flow 
downstream  of  an  abrupt  1:2  expansion .  Measurements  of  mean 
and  RMS  velocities  were  performed  in  a  water  flow  with  a 
laser-Doppler  anemometer.  In  the  upstream  tube,  the 
Reynolds  number  was  varied  from  30,000  to  100,000  and  the 
swirl  number  was  varied  from  zero  to  1.2.  For  low  swirl 
levels,  as  the  core  flow  passed  through  the  expansion,  it 
departed  the  axis  of  symmetry  and  precessed  about  that  axis 
at  frequencies  on  the  order  of  1  Hz.  As  swirl  was  increased 
to  moderate  levels,  the  flow  became  axisymmetric  with  on- 
axis  recirculation  marking  the  onset  of  vortex  breakdown. 
At  the  highest  swirl  levels,  flow  on  the  tube  centerline  was 
in  the  same  direction  as  the  mean  flow,  with  reverse  flow 
occurring  Just  off-axis .  Turbulence  intensities  at  the 
highest  swirl  levels  were  found  to  reach  60k.  As  the  swirl 
was  increased  from  zero  to  its  maximum  value,  the  flow 


reattachment  point  moved  upstream  from  9  to  2  step  heights. 


NOMENCLATURE 


Diaaeter  of  upstreaa  tube 
Oiaaeter  of  downstreaa  tube 


Precession  frequency  of  PVC 
Step  height;  (D2  -  D)/2 
Turbulence  kinetic  energy 
Precessing  Vortex  Core 


Voluaetric  flow  rate 


Radial  coordinate 


Radius  of  upstreaa  tube 
Radius  of  downstreaa  tube 


Reynolds  nuaber  in  upstreaa  tube;  UD/v 

Swirl  nuaber  in  upstreaa  tube,  see  equation  (1) 

RMS  velocity  noraalized  with  Ux  (*  100) 

Local  Bean  axial  velocity 


Axial  velocity  averaged  over  cross  section 

Maxiaua  axial  velocity  in  upstreaa  tube 

Local  aean  tangential  velocity 

Axial  distance  froa  expansion  face 

Reattachaent  length 

Thickness  of  viscous  sublayer 

Ratio  of  D/D2 

Molecular  viscosity  of  fluid 
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INTRODUCTION 


Turbulent  swirling  flow  through  an  abrupt  axisyaaetric 
expansion  is  a  coaplex  flow  possessing  several  distinctly 
different  flow  regiaes,  either  one  or  two  recirculation 
regions,  extreaely  high  levels  of  turbulence,  and  periodic 
asyaaetries  under  soae  conditions.  An  accoapanying 
elevation  of  heat  transfer  rates  is  a  principle  aotivation 
for  the  addition  of  swirl  to  flows  in  duap  coabustors  of  gas 
turbine  engines  and  in  solid  fuel  raajet  coabustors.  The 
objective  of  the  present  investigation  was  to  exper iaental ly 
exaaine  these  flowfields  in  soae  detail. 

The  sudden-expansion  geoaetry  produces  aixing  rates 
downstreaa  of  the  expansion  which  are  substantially  higher 
than  those  which  would  be  obtained  at  the  sane  Reynolds 
nuaber  in  the  entrance  region  of  a  pipe.  This  enhancement 
in  nixing  occurs  in  spite  of  a  recirculation  region 
extending  about  nine  step  heights  downstreaa  froa  the 
expansion.  In  this  recirculation  region.  Bean  velocities 
are  typically  only  ten  percent  as  high  as  those  found  in  the 
core  flow.  The  elevated  aixing  rates  are  due  to  very  high 
levels  of  turbulence  kinetic  energy  generated  by  shearing  as 
the  core  flow  issues  into  the  larger  pipe.  Near  the  tube 
wall,  where  length  scales  are  saall,  dissipation  doainates 
since  the  dissipation  is  inversely  proportional  to  the 
length  scale.  But  in  the  high-shear  regions  away  froa  the 
wall,  length  scales  are  large  and  dissipation  rates 
consequently  low.  Thus,  turbulence  kinetic  energy  generated 


in  the  shear  layer  dissipates  relatively  slowly  and  its 
levels  are  auch  higher  than  would  be  found  in  ordinary  pipe 
flow  where  no  such  internal  shear  layer  exists.  High 
turbulence  kinetic  energy  levels  also  cause  the  thickness 
of  the  (aolecular)  viscosi ty-doainated  sublayer  to  be 
reduced.  Specifically,  for  flows  where  the  principle  energy 
generation  is  not  at  the  wall,  but  rather  reaoved  froa  it  as 
in  the  sudden-expansion  flowfield,  Spalding1  suggested  that 
the  viscous  sublayer  thickness  (y)  changes  with  the 
turbulence  kinetic  energy  (k)  so  that  the  sublayer  Reynolds 
nuaber  (ykl/2/v)  is  a  universal  constant.  Hence,  the 
sublayer  will  becoae  thinner  with  increasing  levels  of 
turbulence  kinetic  energy. 

There  has  been  speculation  (Johnston2)  that  a  snail, 
counter-rotating  corner  eddy  lies  very  close  to  the  face  of 
the  expansion.  Mean  velocities  in  the  corner  eddy  are  of 
the  order  0.01U1  according  to  Johnston,  but  there  are 
apparently  no  velocity  neasureaents  or  flow  visualization 
results  in  the  literature  which  either  confirn  this  value, 
or  support  the  existence  of  this  feature  in  the 
axisyaaetric-expansion  flowfield.  However,  in  a  series  of 
heat  transfer  neasureaents,  Baughn,  Hoffaan,  Takahashi,  and 
Launder3  speculated  that  snail  and  consistent  ainiaa  in 
Nusselt  nuabers  near  the  face  of  the  expansion  were  possible 
evidence  for  the  presence  of  a  corner  eddy.  The 


hypothesized  corner  eddy  is  likely  to  continue  to  defy 
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instruaentation  consists  of  coapsratively  large  probes  or 
probe  voluaes  for  detection  of  such  a  saall  feature  in  this 
restrictive  geoaetry. 

Adding  swirl  to  the  sudden-expansion  flowfield  causes 
an  increase  in  the  width,  growth  rate,  entrainaent,  and 
decay  of  the  core  flow  eaanating  froa  the  upstreaa  tube.  It 
is  also  found  that  on-axis  recirculation  (known  as  vortex 
breakdown4)  aay  occur  for  sufficiently  high  swirl  strengths. 
This  recirculation  is  driven  by  an  adverse  pressure-gradient 
on  the  tube  centerline  which  results  froa  the  viscous 
dissipation  of  the  tangential  velocity  coaponent  as  the  flow 
proceeds  downstreaa.  As  swirl  strength  is  increased  froa 
zero,  the  vortex  breakdown  aay  first  be  seen  as  an  on-axis 
ellipsoid  of  recirculating  fluid.  As  the  degree  of  swirl  is 
further  increased,  the  ellipsoid  aay  stretch  in  the 
downstreaa  direction  and  fora  a  tube  of  recirculating  fluid, 
at  least  in  the  sudden-expansion  geoaetry. 

A  further  coaplex  and  little-understood  phenoaenon 
which  frequently  occurs  in  swirling  flows  is  the  existence 
of  an  unsteady  (although  usually  periodic)  asyaaetry  in  the 
flowfield.  These  asyaaetries  are  usually  associated  with 
the  vortex  breakdown  phenoaenon  and  on-axis  recirculation4 . 
Consequently,  they  are  usually  observed  at  aoderate-to-large 
swirl  strengths.  However,  the  present  investigation 
docuaents  an  asyaaetry  like  that  observed  by  Hallett  and 
Gunther*,  which  occurs  at  low  swirl  strengths  in  the  absence 
of  on-axis  recirculation.  This  latter  flow  asyaaetry  is 


characterized  by  the  vortex  emanating  froa  the  upstreaa  tube 
departing  the  axia  of  syaaetry  and  then  preceaaing  about 
that  axis.  This  feature  will  be  referred  to  here  as  the 
precessing  vortex  core  (PVC)  after  Gupta,  et  al* . 

Analytical  prediction  of  the  present  flowfield  is 
sufficiently  coaplex  that  it  is  aanageable  only  if  the  flow 
is  assuaed  to  be  steady  and  axisyaaetric.  With  these 
siaplif ications  there  is  no  potential  for  predicting  the 
unsteady  three-diaensional  asyaaetry  which  occurs. 
Furtheraore,  Sultanian's7  recent  coaputations  of  this  flow 
had  difficulty  in  accurately  predicting  the  extent  of  the 
on-axis  recirculation  zone  and  turbulence  intensities 
downstreaa  of  the  expansion.  Sultanian  also  found  his  aodel 
to  be  quite  sensitive  to  the  inlet  conditions,  especially 
turbulence  intensity.  With  this  in  Bind,  we  present 
aeasureaents  upstreaa  of  the  sudden  expansion  to  facilitate 
subsequent  aodeling  efforts. 

At  this  point,  it  is  convenient  to  define  several 
scales  and  independent  variables  which  will  be  used  in  the 
following  discussion.  There  are  two  length  scales  required 
in  the  axisyaaetric  sudden-expansion  problea.  The  first  is 
the  step  height  (h)  which  experience  has  shown  to  be 
reasonably  well  suited  for  correlation  of  reattachaent 
lengths.  A  second  necessary  length  scale  is  either  the 
upstreaa  or  the  downstreaa  tube  diaaeter.  Here,  the 
upstreaa  tube  diaaeter  (D)  is  eaployed.  The  Reynolds  nuaber 
is  based  on  the  diaaeter  of  the  upstreaa  tube  and  the 
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average  velocity  in  the  upstreaa  tube.  Swirling  flows  are 
coaaonly  characterized  by  the  following  definition  for  a 


device-independent  swirl  nuaber: 
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The  swirl  nuaber  aay  be  physically  interpreted  as  the  ratio 
of  axial  fluxes  of  swirl  and  linear  aoaentun,  divided  by  a 
characteristic  radius. 


PREVIOUS  INVESTIGATIONS 

Sudden  Expansion  Flow  without  Swirl.  Axial  flow 
through  a  sudden  axisyaaetric  expansion  is  a  fairly  well 
studied  problea8-13  which  represents  the  liaiting  case  of 
zero  swirl  against  which  current  results  can  be  coapared. 
The  widely  referenced  set  of  data  by  Chaturvedi8  includes 
aean  velocities  and  turbulence  quantities  aeasured  with  a 
hot  wire  aneaoaeter.  However,  a  check  of  Chaturvedi’s  aass 
balances  yields  prof i le-to-prof i le  variations  as  high  as 
30*.  Chaturvedi  atteapted  to  saooth  the  data,  but  the 
arbitrary  nature  of  this  correction  reduces  one's  confidence 
in  it.  Freeaan*  used  a  laser-Doppler  aneaoaeter  (LOA)  to 
aeasure  axial  aean  velocity  and  turbulence  intensity  while 
Moon  and  Rudinger10  report  only  an  axial  aean  velocity  froa 
their  LOA  aeasureaents .  Yang  and  Yu11  report  turbulence 
quantities  and  aean  velocities,  also  obtained  with  an  LOA, 
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but  the  validity  of  tbeir  data  has  been  called  into  question 
recently14  due  to  significant  aass  balance  discrepancies 
which  are  actually  higher  than  those  quoted  in  the  paper. 
Among  the  various  studies,  the  measurements  of  Gould, 
Stevenson  A  Thompson12  and  Stevenson,  Thompson  A  Gould13 
appear  to  be  the  aost  complete  and  most  closely  related  to 
the  present  work. 

The  present  state  of  the  art  of  computational  flow 
modeling  is  such  that  the  sudden-expansion  problem  (purely 
axial  flow)  is  now  fairly  well  handled  by  various  schemes. 
The  reader  is  referred  to  Gossan,  Khalil,  and  Whitelaw19  and 
Stevenson,  Thompson,  and  Gould13  for  discussions  of  k- « 
modeling,  and  to  Minh  and  Chassaing1*  and  Sultanian7  for 
application  of  Reynolds  stress  modeling  to  this  problem. 

Axisymmetric  Expansion  with  Swirl.  There  have  been 
several  recent  investigations  reporting  measurements  in 
swirled  flows  through  sudden  expansions17*20,  but  all  have 
used  intrusive  probes,  even  though  it  is  known  that  such 
probes  can  significantly  alter  flowfields  with 
recirculation.  In  fact,  these  studies  are  concerned 
primarily  with  the  development  of  measurement  techniques 
using  five-hole  pitot  tubes  and  hot  wire  anemoaetry  for 
application  in  multi-dimensional  complex  flows. 
Consequently,  these  four  papers  might  be  considered  work  in 
progress  on  the  development  of  measurement  techniques  rather 
than  a  collection  of  results  available  for  comparison 


purposes . 


Vortex  Breakdown  and  the  PVC.  There  have  been  a 
nuaber  of  analytical  investigations  of  vortex  breakdown  (see 
reviews  by  Hall21  and  Leibovich22)  but  the  asynaetries  in 
swirled  flows  are  so  coaplex  and  irregular  that  these 
analytical  treataents  have  been  aostly  unsuccessful.  Thus 
the  priaary  body  of  inforaation  about  unsteady  asyaaetries 
in  swirling  flows  has  been  gathered  in  experiaental 
studies4 * 23-26  .  The  flow  geoaetries  in  these  experiaents 
are  all  axisyaaetric ,  but  a  wide  variety  of  configurations 
are  represented.  These  include  straight  tubes23 , 
diffusers4 ■ 25 • 26 ,  sudden  contractions23,  and  unconfined 
swirling  jets24.  Although  the  geoaetries  are  diverse,  the 
nature  of  the  asyaaetrical  flows  observed  in  the  various 
experiaents  is  reaarkably  siailar.  The  single  feature 
coaaon  to  all  of  these  flows  is  a  precession  of  the  flow 
about  the  tube  axis  in  conjunction  with  vortex  breakdown. 

All  of  the  prior  vortex  breakdown  experiaents  were 
essentially  flow  visualization  studies.  In  those  cases 
where  hot  wire  and  laser-Doppler  aneaoaetry  were 
eaployed4 • 2 3 • 2 4 ,  they  were  used  to  look  for  a  sinusoidal 
variation  in  aean  velocity  as  the  asyaaetry  swept  past  the 
point  of  aeasureaent.  Cassidy  and  Falvey23  noted  that  the 
precession  frequency  was  independent  of  Reynolds  nunber  for 
Re  &  10*.  In  very  careful  flow  visualization  studies,  Faler 
and  Leibovich4  identified  six  distinct  disturbance  nodes 
whose  flow  regiaes  could  be  characterized  by  Reynolds  and 


thea  to  conclude  that  there  are  ao  truly  axisyaaetric 
disturbaace  patterns  ia  these  flows. 

There  is  aa  important  distinction  between  the  papers 
aentioned  above  and  the  work  of  Hallett  and  Gunther5.  In 
the  previously  aentioned  studies,  the  flow  asyanetries  occur 
only  in  conjunction  with  vortex  breakdown,  but  Hallett  and 
Gunther’s  PVC  in  a  sudden  expansion  occurs  only  for  swirl 
strengths  below  those  associated  with  vortex  breakdown.  In 
fact,  with  increasing  swirl,  the  periodicity  of  the  PVC 
becaae  weaker  and  less  distinct  until  just  before  onset  of 
vortex  breakdown,  it  disappeared  altogether.  Also,  they 
noted  that  precession  was  strongest  and  aost  regular  at  low 
swirl,  while  at  higher  swirl  the  aotion  becaae  increasingly 
irregular.  Hallett  and  Gunther  observed  that  the  aaplitude 
of  the  PVC  dissipated  with  increasing  downstreaa  distance  as 
the  PVC  becaae  aore  coincident  with  the  tube  axis.  Finally, 
they  report  no  evidence  of  flow  asyaaetry  upstreaa  of  the 
expansion  in  velocity  aeasureaents  Bade  with  a  five-hole 
pitot  probe. 

EXPERIMENTAL  APPARATUS 

A  water  flow  loop,  constructed  of  stainless  steel,  and 
shown  scheaat ical ly  in  Figure  1,  coaprised  the  sain  element 
of  the  test  facility.  Swirl  was  generated  by  supplying  a 
variable  portion  of  the  flow  through  tangential  slots  as 
indicated  by  Figure  2.  Inside  diaaeters  of  the  axial  inlet 


tube,  the  swirler  insert,  and  the  upstream  test  section  were 
5.08  ca.  The  axial  inlet  tube  was  31  dianeters  long  to 
allow  axial  flow  developaent,  and  the  sudden  expansion  was 
15  diaaeters  downstreaa  of  the  swirl  generator.  Flow  rates 
to  the  slots  and  the  axial  inlet  tube  could  be  controlled 
independently,  thus  providing  the  capability  to  continuously 
vary  swirl  strength.  Flow  rates  were  Measured  with  turbine- 
type  flowaeters.  An  in-line  filter  was  used  to  reaove 
particles  noainally  larger  than  1  froa  the  water. 

The  tube  upstreaa  of  the  expansion  was  Bade  of 
Plexiglas  which  was  bored,  honed,  and  polished  to  a  final  ID 
of  5.078  ±.008  ca.  To  allow  LDA  Measurements  close  to  the 
expansion,  the  tube  and  attached  expansion  face  extended 
into  the  downstreaa  tube  so  that  structural  flanges  and 
bolts  did  not  interfere  with  the  laser  beams.  Measurements 
were  thus  possible  1  ca  downstreaa  of  the  expansion.  The 
downstreaa  tube  was  not  Machined  or  honed  due  to 
complications  associated  with  its  relatively  large  size. 
Consequently,  it  was  very  slightly  oval  with  an  ID  of  9.850 


±.020  ca  and  00  of  10.767  ±.003  ci 


Thus  the  expansion 


ratio  was  1:1.94.  The  length  of  the  downstreaa  tube  was 
1.04  a. 


The  LDA  Optical  System.  The  laser-Doppler  anemometer 
was  a  conventional,  single-component  system  operated  in  the 


dual-beaa  Mode.  The  systea  included  a  15  aW  He-Ne  laser,  a 
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components  produced  an  ellipsoidal  probe  volume  whose 
nominal  1/e2  extent  was  .91  mm  long  and  .09  mm  in  diameter. 

Both  transmitting  and  forward-scatter  receiving  optics 
were  mounted  on  a  single  aluminum  channel  which  in  turn  was 
rigidly  affixed  to  the  table  of  a  three-axis  milling 
machine.  To  obtain  the  desired  1  a  travel  in  the  axial 
direction,  the  milling  table  was  permanently  fastened  to 
four  precision  linear  bearings  which  rode  on  two  parallel 
steel  shafts.  A  dial  indicator  was  used  to  monitor  the 
radial  position  of  the  probe  volume. 

EXPERIMENTAL  PHOCBDURBS 

Artificial  seeding  of  the  flow  was  not  required.  The 
test  loop  was  filled  with  tap  water  having  a  relatively  high 
mineral  and  particulate  content.  The  water  was  filtered 
briefly  after  each  filling  of  the  loop  to  remove  particles 
nominally  larger  than  1  p m.  Filtering  to  this  size  was 
consistent  with  the  Melling  and  Whitelaw27  suggestion  that 
particulates  smaller  than  10  pm  will  adequately  follow  the 
flow  up  to  frequencies  of  500  Hz.  For  most  water  flows,  the 
bulk  of  the  energy-containing  eddies  have  frequencies  in 
this  range. 

Both  axial  and  tangential  components  of  mean  velocity 
and  RMS  turbulence  levels  were  measured  on  a  dense  grid  of 
points  lying  in  a  horizontal  plane  through  the  tube 
centerline.  Included  in  the  grid  of  measurement  stations 


were  two  upstreaa  locations  at  X/D  =  -2.0  and  -0.5. 
Locations  for  profiles  in  the  downstream  tube  were  chosen  to 
optimize  resolution  in  the  near-expansion  region  where 
velocity  and  turbulence  levels  change  rapidly  with  X/D. 
Corrections  for  optical  refraction  of  the  laser  beaas  at  the 
air-Plexiglas  and  Plexiglas-water  interfaces28  were  employed 
to  locate  the  probe  voluae  at  even  intervals  in  the  radial 
direction.  For  aost  cases  with  swirl,  profiles  were  aade 
across  the  entire  tube  to  check  for  flow  asyaaetry  even 
though  asyaaetry  was  only  found  in  the  subcritical-swirl 
flows . 

Measurements  at  radius  ratios  ±0.95  were  attempted  for 
all  profiles.  However,  measurement  close  to  walls  is 
generally  difficult  since  scattered  light  froa  the  walls 
results  in  poor  signal-to-noise  ratios.  In  the  present 
work,  spurious  wall  reflections  were  largely  overcome  by 
collecting  scattered  light  at  about  5  degrees  off  the 
forward  scatter  axis,  thus  truncating  the  probe  voluae 
slightly  and  keeping  the  light-collection  optics  out  of  the 
horizontal  plane  which  contains  aost  of  the  disruptive  stray 
light.  Using  this  technique,  results  at  r/Ra  =  ±.95  in  the 
downstream  tube  were  consistent  and  credible.  In  the 
upstreaa  tube,  because  of  its  saaller  size,  the  saae  degree 
of  credibility  extended  only  to  r/R  =  ±.90. 

Velocity  biasing  was  eliminated  by  randoa  saapling  in 
the  present  experiaents.  Durao,  Laker,  and  Whitelaw3* 
suggest  that,  for  data  collection  percentages  of  less  than 


40*,  the  average  velocity  obtained  will  be  less  than  2* 


higher  than  the  true  aean  velocity.  Stevenson,  et  el.13 
suggest  that  the  velocity  bias  will  be  effectively 
eliminated  for  collection  percentages  on  the  order  of  1*. 
For  this  work,  a  computer  sampled  the  output  from  a  counter 
processing  device  at  a  rate  of  130  Hz.  Data  rates  often 
fell  to  about  4000  Hz  at  the  near-wall  grid  points  of 
r/R2  =  ±.95,  but  usually  ranged  from  8000  to  40,000  Hz 
elsewhere.  Hence  the  worst-case  collection  percentages  were 
about  3*  near  the  tube  walls.  The  waiting  period  also 
minimizes  the  potential  bias  caused  by  a  single  particle 
generating  multiple  measurements  before  leaving  the  probe 
volume. 

Mean  and  RMS  velocities  were  determined  from  sample 
sizes  of  4000  data  points.  The  statistical  error30 
associated  with  this  sample  size  is  ±1*  in  the  mean-velocity 
measurement  for  a  local  turbulence  intensity  (TI)  of  70*  and 
about  ±2*  in  the  measurement  of  TI.  A  worst-case 
computation  of  the  spatial  velocity  biasing  due  to  the 
finite  probe-volume  size31  suggests  that  the  spatially 
averaged  TI  is  only  0.6*  higher  than  the  TI  at  the  probe 
volume’s  center. 

Investigation  of  the  PVC.  The  experimental 
examination  of  the  precessing  vortex  core  and  the  vortex- 
breakdown  bubble  consisted  of  both  flow  visualization  and 
selective  probing  with  the  LDA.  Some  modest  success  in 
visualizing  the  flow  was  obtained  using  air  bubbles  and 


high-intensity  lighting  in  two  different  procedures.  The 
first  and  lore  useful  method  was  to  introduce  approximately 
one  liter  of  air  into  the  250  liter  capacity  of  the  test 
loop.  The  air  and  water  were  then  mixed  by  operating  the 
loop  for  a  short  time.  The  ensuing  air  bubbles  were  so 
small  as  to  be  almost  invisible  to  the  eye  with  ordinary 
room  lighting,  but  with  the  use  of  a  high  intensity 
photographic  light  source,  a  "mist"  of  bubbles  could  be  seen 
well  enough  to  perceive  qualitative  details  of  the  flow. 
These  bubbles  were  sufficiently  small  that  they  showed  no 
discernible  tendency  to  either  rise  to  the  top  of  the  tube 
or  collect  on  the  tube  centerline,  but  rather  seemed  to 
follow  the  flow.  When  the  flow  field  was  visualized  in  this 
fashion  it  was  often  difficult  to  determine  what  was 
happening  in  the  tube.  In  particular,  while  it  was  clear 
that  the  vortex  from  the  upstream  tube  was  entering  the 
downstream  tube  asymmetrically  and  precessing,  vigorous 
activity  in  the  near-wall  recirculation  zone  complicated  the 
examination  so  that  the  direction  of  precession  at  very  low 
swirl  numbers  (S  1  .15)  could  not  be  determined.  To  aid  in 
the  resolution  of  this  dilemma,  air  was  injected  through  a 
small  total-pressure  probe  on  the  centerline  of  the  upstream 
tube.  The  air-injecting  probe  was  located  just  downstream 
of  the  swirl  generator.  As  the  larger  bubbles  which  were 
produced  in  this  way  moved  downstream,  they  were  pinned  on 
the  tube  centerline  by  centrifugal  forces.  As  they  passed 


through  the  expansion,  they  Barked  the  vortex  axis  and  thus 
revealed  the  direction  of  vortex  precession. 

Precession  frequency  inforaation  was  gathered  by 
aonitoring  the  counter-processor's  analog  output  on  both  a 
strip-chart  recorder  and  a  spectrua  analyzer  to  obtain  a 
real-tine  variation  of  aean  velocity.  An  RMS  voltaeter  with 
adjustable  tiae  constant  was  connected  between  the  counter's 
output  and  the  recorder  or  analyzer  so  that  the  higher 
frequencies  associated  with  turbulent  fluctuations  could  be 
filtered  out. 

The  nethods  of  Kline  and  McClintock32  were  eaployed  to 
deternine  that  the  largest  uncertainties  were  about  2X  in 
Reynolds  nuaber,  8 X  in  swirl  nuaber,  10*  in  Strouhal  nuaber, 
and  1*  in  probe  volune  positioning.  Uncertainties  in  aean 
and  RMS  velocities  due  to  the  aany  possible  biases  and 
broadening  errors  are  estimated  to  be  about  ±3*  and  tlO*, 
respectively. 

RESULTS  AND  DISCUSSION 

Velocity  and  TI  Distributions.  The  results  for  U,  u’  , 
and  v'  for  unswirled  flow  at  Reynolds  nuabers  of  30,000, 
60,000,  and  100,000  are  shown  in  Figure  3.  The  aean  and 
fluctuating  velocities  in  Figure  3  have  been  noraalized  with 
the  axial  centerline  velocity  occurring  in  the  upstream 
tube.  Mean  velocities  for  the  three  Reynolds  nuabers 
collapse  to  single  curves,  but  when  exaaining  the  regions 


far  downstream  (X/D  =  18)  it  can  be  seen  that  the  turbulent 
fluctuations  have  apparently  dissipated  faster  for  larger 
Reynolds  numbers  (the  TI  for  Re  =  60,000  and  Re  =  100,000 
were  virtually  identical).  Figure  3  indicates  that  the 
axial  TI  has  decreased  to  a  nominal  background  level  of 
about  2.5%  for  Re  =  100,000  while  it  remains  near  8%  for 
Re  =  30,000.  The  difference  in  TI's  is  due  to  higher  rates 
of  dissipation  at  the  larger  Reynolds  numbers.  Dissipation 
usually  scales  as  u’3/J,  or  ( u ’ /U) 3 / ( J/U3 ) ,  so  that  for 
length  scales  (J)  and  turbulence  intensities  (u’/U  a  100)  of 
the  sane  order,  dissipation  increases  with  increasing  U,  or 
increasing  Reynolds  number. 

For  unswirled  flows,  we  see  that  a  state  of  near¬ 
isotropy  in  u’  and  v’  exists  at  X/D  -  -2.0,  and  then  again 
far  downstreaa  after  the  flow  has  redeveloped.  However, 
throughout  auch  of  the  interaediate  region  the  axial  TI  is 
approximately  30%  greater  than  the  tangential  TI.  The  peak 
values  for  both  of  these  quantities  are  generally  coincident 
and  lie  in  the  region  bounded  by  the  edge  of  the  shear  layer 
and  the  tube  centerline.  For  each  Reynolds  number,  these 
maximum  values  are  on  the  order  of  20  -  22%  for  axial  TI  and 
about  15%  for  tangential  TI.  This  aaxiaua  value  (and  the 
distributions  of  U  and  u')  compare  very  favorably  with  the 
work  of  Stevenson,  Thompson,  and  Gould13  who  report  a 
aaxiaua  axial  TI  of  22%. 

Swirled-flow  data  for  five  supercr i t ical -swirl  cases 
are  shown  in  Figures  4-8.  For  these  highly  swirled  flows, 


the  aaxiaua  axial  velocity  in  the  upstreaa  tube  (U1) 
occurred  near  r/R  *  .8,  and  it  is  this  value  which  is  used 
for  the  noraalization  of  aean  and  RMS  velocities.  Figures 
4-8  deaonstrate  a  large  influence  of  the  downstreaa  flow  on 
the  X/D  =  -0.5  profiles  for  all  supercritical-swirl  cases. 
The  influence  is  especially  strong  on  the  TI’s  whose 
aagnitudes  and  distributions  are  changed  draaatically  from 
the  X/D  =  -2.0  station.  The  turbulence  intensities  continue 
to  be  highly  non-isotropic  in  the  downstreaa  tube  for  all 
the  supercritical-swirl  flows.  These  results  have  iaportant 
raaif ications  on  aodeling  of  this  flow  since  they  iaply  the 
frequently  used  k-s  aodel,  with  its  assuaption  of  isotropy, 
will  be  unable  to  accurately  predict  the  observed  features 
of  highly  swirled  flows. 

The  highest  swirl  strength  achieved  in  this  set  of 
experiaents  is  shown  in  Figure  8,  where  we  see  two  features 
not  present  in  the  flows  with  lower  swirl  nuabers.  The 
first  is  that,  as  the  flow  developaent  proceeds  downstreaa 
froa  X/D  =  4,  the  centerline  velocity  is  positive  with 
recirculation  just  off-axis.  We  also  find  that  the  on-axis 
tangential  velocity  gradient  is  steepening  throughout  the 
downstreaa  tube,  consequently  producing  greeter  shearing  and 
ever-increasing  turbulence  intensities.  Both  of  these 


trends  continued  through  X/D  =  18, 


These  two  features 


cannot  be  wholly  discerned  froa  the  data  presented  for 
Re  -  60,000,  S  =  1.16,  but  froa  Figure  6  and  aeasureaents 
aade  further  downstreaa  (those  for  X/D  >  10  not  shown  here), 
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the  trends  toward  positive  centerline  velocity  and 
continuously  increasing  TI  do  sees  to  be  present  and  hence 
consistent  with  the  Re  =  100,000,  S  =  1.23  data.  That 
centerline  TI’s  continue  to  increase  to  the  end  of  the  test 
section  is  a  limitation  of  the  test  section's  length.  In  a 
longer  tube  the  decay  of  swirl  would  give  rise  to  a  maximum 
in  TI  at  soae  axial  location.  However  it  is  perhaps 
surprising  that  this  flow  condition  with  its  high  diffusion 
rates  and  relatively  short  wall  recirculation  zone  is  still 
evolving  so  far  downstream  at  X/D  =  18. 

For  the  swirled  flows  in  general,  the  peak  value  for 
axial  TI  always  occurs  in  the  shear  layer  near  r/R2  =  0.5, 
and  the  maximum  value  of  tangential  TI  is  always  found  along 
the  tube  centerline.  For  the  swirled  flows  there  is  a 
considerable  divergence  in  the  behaviors  of  axial  and 
tangential  turbulence  intensities  for  X/D  =  -0.5  and 
throughout  the  downstream  tube.  It  can  be  seen  from 
Figures  4-8  that,  along  the  tube  centerline,  the  tangential 
TI  is  typically  twice  the  axial  TI.  At  the  same  time  in  the 
shear  layer  around  r/R2  =  0.5,  the  axial  TI  is  commonly 
twice  the  tangential  TI.  The  largest  axial  TI  are  on  the 
order  of  45%  for  moderate  swirl  (.60  <  S  <  .77)  and  58  -  65% 
for  high  swirl  (.98  <  S  <  1.23).  In  each  case  the 
corresponding  maximum  tangential  TI  is  always  several 
percent  less  than  the  axial  value.  For  the  swirled  flows 
these  msxiaum  values  were  found  near  X/D  =  0.5,  while  for 
the  unswirled  flows,  they  were  found  between  3  <  X/D  <  6. 


A  composite  of  axial  centerline  velocities  is  shown  in 


Figure  9.  For  the  unswirled  flows,  the  flat  portion  of  the 


curve  for  X/D  >  12  indicates  that  the  velocity  profile  has 


redeveloped. 


From  an  area-ratio  argument,  one  would 


anticipate  that  for  the  present  expansion  ratio  of  0.51  that 


the  normalized  centerline  velocity  would  achieve  a 


downstreaa  value  of  (0.51)2  or  0.265.  The  actual  value 


reached  was  a  consistent  0.24,  which  suggests  that 


downstreaa  profiles  are  flatter  and  aore  developed  than 


those  upstream  of  the  expansion.  That  this  is  true  can  be 


seen  in  the  velocity  profile  plots  of  Figure  3.  It  is 


hypothesized  that  the  profile  upstream  of  the  expansion  is 


not  quite  fully  developed  due  to  slight  tube-wall 


irregularities  associated  with  the  swirl  generator’s  slots, 


pipe  joints  located  at  the  swirl  generator,  and  a  flanged 


pipe  connection  approximately  four  diaaeters  upstream  of  the 


X/D  =  -2.0  measuring  station. 


The  normalized  centerline  velocity  has  no  Reynolds- 


number  dependence  for  the  unswirled  flows.  If  we  assume 


this  also  to  be  true  for  swirled  flows,  then  Figure  9 


suggests  that  there  is  a  swirl  number  between  0.74  and  0.98 


which  gives  a  maximum  reverse  velocity. 


Further,  the 


laximum  reverse  velocity,  which  is  seen  to  occur  between 


X/D  -  0.5  and  X/D  =  1.3,  decreases  in  magnitude  as  swirl 


nuaber  is  increased  beyond  S  =  .98.  This  suggests  that 


there  aay  be  swirl  nuabers  greater  than  1.23  for  which  the 


centerline  velocity  aay  always  be  positive,  as  for  the 
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unswirled  flows.  Unfortunately,  S  =  1.23  was  the  upper 
limit  on  swirl  strength  available  from  the  equipment  used  in 
this  study,  so  that  this  hypothesis  could  not  be 
investigated  further. 

For  the  present  work,  the  position  of  U  =  0  points  at 
radial  locations  of  r/R2  =  .8,  .9,  and  .95,  was  determined 
by  linear  interpolation  between  adjacent  grid  points  for 
which  U  had  opposite  signs.  The  U  =  0  points  froa  the  three 
radial  locations  were  then  fitted  with  a  spline  and  the 
resulting  curve  extrapolated  to  the  wall  to  find  the 
reattachaent  "point".  Reattachaent  lengths  obtained  in  this 
fashion  for  the  present  data  are  shown  in  Table  1.  For  the 
three  unswirled  cases,  the  reattachaent  lengths  agree  well 
with  those  reported  for  unswirled  flows  in  the  previously 
mentioned  investigations  (also  shown  in  Table  1).  We  know 
of  no  existing  data  which  can  be  used  for  comparison  of 
reattachaent  lengths  in  swirled  flows.  When  reattachaent 
lengths  are  plotted  against  swirl  nuaber,  the  resulting 
curve  is  independent  of  Reynolds  number  and  appears  to  be 
asyaptotic  to  xr/h  *  1.7  as  swirl  number  increases  beyond 
1.2.  The  reattachaent  lengths  given  here  are  actually 
average  values  determined  from  four  to  eight  sets  of 
velocity  data  obtained  in  the  region  of  reattachaent.  The 
slight  variations  found  in  these  velocities  were  sufficient 
to  appreciably  alter  the  computed  lengths,  this  being 
especially  true  for  the  unswirled  flows  where  variations  in 
reattachaent  length  of  t0.5  step  heights  were  not  uncommon 


between  individual  data  seta.  The  variation  in  reattachaent 
length  for  the  swirled  cases  was  a  aore  aodest  tO.l  step 
heights  for  individual  data  trials.  It  is  probably 
reasonable  to  consider  these  variations  as  representative  of 
randoa  fluctuations  in  the  reattachaent  zone’s  width  rather 
than  as  an  uncertainty  in  the  aeasureaent,  but  it  is 
difficult  to  separate  these  two  effects. 

Mass  balances  obtained  from  integration  of  velocity 
profiles  have  becoae  a  coaaonly  used  standard  for  appraising 
the  credibility  of  internal-flow  velocity  data.  For  a 
particular  flow  condition  in  the  present  work,  the  largest 
differences  between  the  aass  flux  at  any  one  profile  and  the 
average  aass  flux  for  all  profiles  (at  that  flow  condition) 
were  between  3*  and  5*.  The  locations  of  poorest  agreeaent 
were  randoaly  scattered  in  the  axial  direction.  Tabular 
velocity  data  and  further  details  are  available  in  the 
dissertation  by  Del lenback3 3 . 

The  PVC  and  Vortex  Breakdown.  Information  which  can 
be  generalized  about  swirling-flow  asyaaetries  and  the  PVC 
froa  previous  studies  is  very  sparse,  largely  due  to  the 
coaplexity  of  the  flow’s  structure.  Historically,  the  aost 
easily  and  coaaonly  aeasured  feature  of  these  asyaaetries 
has  been  the  relationship  between  the  swirl  nuaber  and 
precession  frequency.  This  relationship  also  constitutes 
the  principle  result  of  the  present  study  of  the  PVC.  While 
the  coaputation  of  swirl  nuaber  requires  knowledge  of  the 
aean  velocity  profiles,  which  are  not  known  a  priori,  the 


•wirl  number  can  also  ba  related  to  the  ratio  of  lass  fluxes 
entering  the  swirl  generator.  For  weakly  swirled  flows 
(S  <  .15),  the  swirl  ousber  was  obtained  fron  an  algebraic 
relationship6  between  the  axial  and  tangential  sass  fluxes 
which  results  fron  assuming  plug  flow  with  superimposed 
solid-body  rotation.  A  second,  experimentally  determined 
relationship  between  the  mass-flow  ratio  and  the  swirl 
number  (from  integrated  velocity  profiles)  was  deduced  for 
higher  swirl  numbers.  Precession  frequencies  were  examined 
only  at  the  two  limiting  Reynolds  numbers  of  30,000  and 
100,000. 

Precession-frequency  data  are  combined  with  flow- 
visualization  observationa  to  summarize  vortex  breakdown  and 
the  PVC’s  swirl-dependent  behavior  in  Table  2.  Previous 
investigations  have  found  that  swirling  flow  asymmetries 
usually  precess  in  the  same  direction  as  the  mean  swirl  for 
confined  flows,  and  in  the  opposite  direction  for  free  jets. 
However,  Escudier25  reports  that  his  sudden  expansion  data 
contradict  this  general  rule.  In  the  present  experiments, 
the  PVC  precessed  with  the  mean  swirl  for  larger  swirl 
numbers,  and  against  the  mean  swirl  for  low  swirl  numbers. 
Table  2  shows  that,  while  results  for  the  two  Reynolds 
numbers  are  similar,  specific  events  occur  at  slightly 
different  swirl  numbers.  Two  regions  were  especially 
difficult  to  resolve.  These  were  the  precession  frequencies 
for  very  low  swirl  (S  <  .1)  and  the  swirl  number  at  which 
the  direction  of  precession  changes.  On  the  other  hand,  two 


points  which  had  very  sharp  transitions  as  the  swirl  nunber 
was  changed  were  the  transition  fro*  PVC  to  vortex 
breakdown,  and  the  transition  froa  a  bubble-like  vortex 
breakdown  to  a  full-length  tube  of  recirculating  fluid  on 
the  tube  centerline.  Throughout  the  regine  of  swirl  numbers 
for  which  the  vortex-breakdown  bubble  exists,  considerable 
unsteadiness  of  the  bubble’s  location  was  noted  in  both  the 
flow  visualization  and  LDA  data.  This  unsteadiness  aade  it 
impossible  to  accurately  measure  either  the  velocities 
inside  the  bubble  or  the  bubble's  dimensions,  but  it  was 
noted  that  the  bubble  extended  to  about  2-1/4  D  downstream 
of  the  expansion.  Downstream  of  the  bubble,  it  was  observed 
that,  while  fluid  on  the  tube  axis  was  mostly  stagnant  in 
the  mean,  it  was  also  quite  unsteady,  sometimes  showing  a 
tendency  to  drift  randomly  either  upstream  or  downstream. 
When  the  PVC  was  present,  flow  oscillations  could  be 
detected  at  the  upstream  station  of  X/D  =  -0.5,  but  not 
further  upstream  at  X/D  =  -2.0.  Although  several  of  the 
investigations  cited  previously  observed  the  PVC  at 
supercritical-swirl  numbers,  the  PVC  in  the  present  study 
could  only  be  detected  at  subcr i t i cal -swi r 1  numbers. 
Finally,  none  of  the  results  shown  in  Table  2  displayed  any 
apparent  hysteresis  in  swirl  number. 

Frequency  information  taken  from  the  stripcharts  (or 
spectrum  analyzer,  which  agreed  closely)  was  converted  to 
Strouhal  numbers23  (fD3/Q,  where  f  is  the  precession 
frequency  and  Q  the  volumetric  flow  rate)  and  used  to 
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generate  Figure  10.  The  scarcity  of  data  points  for  S  <  .1 
and  near  the  zero-frequency  crossovers  is  illustrative  of 
the  previously  aentioned  inability  to  resolve  frequencies  at 
these  transitory  swirl  nunbers.  Away  froa  these  swirl 
nuibers,  saooth  curves  drawn  through  the  available  data 
points  fit  well  and  thus  give  confidence  in  extrapolating 
the  swirl  nuaber  for  the  zero  frequency  crossovers. 
Figure  10  also  shows  clearly  the  dependency  of  the 
precession  frequency  on  Reynolds  nuaber. 

The  frequencies  of  precession  seem  to  follow  a  trend 
consistent  with  other  investigations  of  asymmetries  in 
swirling  flows  as  shown  by  Figures  11  and  12.  Although 
Syred  and  Beer’s24,  and  Cassidy  and  Falvey’s23  experiments 
were  in  dissimilar  geometries  (free  jet;  straight  tube  and 
sudden  contraction,  respectively)  and  the  PVC  occurred  in 
conjunction  with  vortex  breakdown,  there  is  a  noticeable 
congruence  of  results  in  Figure  12.  In  particular,  the 
Strouhal  number  appears  to  be  independent  of  Reynolds  number 
for  Re  >  10s,  but  a  slowly  decreasing  function  of  Reynolds 
nuaber  at  lower  values  of  this  parameter. 


CLOSURE 

The  principle  contribution  of 
been  to  quantify  the  aagnitudes  of 
velocities  as  a  function  of  swirl  and 
abrupt  ax iayaaet r ic-expana ion  flow. 


the  present  study  has 
mean  and  fluctuating 
Reynolds  nuabers  in  an 
In  addition,  a  nuaber 
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of  qualitative  trends  and  features  have  been  observed  which 
■ay  contribute  in  a  general  way  to  the  understanding  of 
complex  shear  flows.  Documentation  of  the  velocity  field  in 
the  upstream  tube  should  be  helpful  to  computational  flow 
modelers,  and  of  general  interest  in  light  of  how  vastly  the 
upstream  flow  changes  just  before  reaching  the  expansion. 
For  low  swirl  levels,  an  unsteady,  three-dimensional  flow 
asymmetry  has  been  observed.  The  asymmetry  is  such  a 
complex  flow  structure  that  only  a  limited  number  of  useful 
measurements  can  be  obtained  which  help  to  specify  its 
structure.  At  the  higher  swirl  levels,  the  flow  becomes 
symmetric  and  extraordinarily  high  levels  of  turbulence  are 
produced.  The  on-axis  recirculation  normally  associated 
with  highly  swirled  flows  changes  such  that  the  on-axis  flow 
is  in  the  downstream  direction  with  recirculation  just  off- 
axis. 
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Figure  12.  Relationship  between  Strouhal  Nueber  and  Reynolds  Nueber 

(Various  Bxperiaents) . 


